Search for CPT Violation with the FOCUS 
Experiment and Measurement of A5 lifetime in the 
decay A5 — > J/i/jA with the D0 Experiment 



Abaz Kryemadhi 



Submitted to the faculty of the Graduate School 
in partial fulfillment of the requirements 
for the degree 
Doctor of Philosophy 
in the Department of Physics, 
Indiana University 



December, 2004 



Accepted by the Graduate Faculty, Indiana University, in partial fulfillment 
of the requirements of the degree of Doctor of Philosophy. 



Doctoral Professor Rick Van Kooten 

Committee (Chairman) 



Professor Alan Kostelecky 



Professor Harold Ogren 



December 7, 2004 Professor Steven Gottlieb 



ii 



Copyright © 2008 
Abaz Kryemadhi 
ALL RIGHTS RESERVED 



iii 



To God: Who created so much beauty in the universe for us to study. 

"The heavens declare the glory of God; the skies proclaim the work of his hands". 

Psalm 19:1 



iv 



Acknowledgments 



First of all, my deepest gratitude goes to my advisors, Dr. Rick Van Kooten and Dr. Rob 
Gardner. Throughout the years of my Ph.D. study, they were always willing to mentor, 
encourage, and support me. Their expertise and enthusiasm in the field of high energy 
physics, as well as their unique insightful way of approaching problems, are among the most 
valuable resources for me. 

I would like to thank Dr. Alan Kostelecky for his explanations of the CPT and Lorentz 
Violation framework and also for his comments and suggestions during the process of pub- 
lishing the CPT paper, which has been used in this thesis. 

I would like to thank Dr. Fred Luehring for his suggestions when I have been at cross- 
roads in my career and his help with some computer software. Along the same lines I would 
like to thank Thorn Sulanke for his help at solving a lot of technical questions. 

I would like to thank people from the FOCUS Collaboration, in particular Jim Wiss, 
Gianluigi Boca, Topher Cowfield, John Link, Eric Vaandering, Kevin Stenson and Harry 
Cheung for their continual help with any question that I had, and also for their input on 
making the first part of this thesis possible. 

I would like to thank Indiana University D0 group, in particular Daria Zieminska and 
Andrzej Zieminski for the usuful disscussions and suggestions during our Indiana University 
D0 group meetings. 

I would like to thank people from the B-Physics group at D0 in particular Vivek Jain, 

v 



Brad Abbott, Rick Jessik, Guenadi Borisov, Andrei Nomerotski, Sergey Burdin for their 
continual help with any question that I had, and also for their input on making the second 
part of this thesis possible. 

I would like to thank friends in the Physics Department, Maciej Swat, Prabudha Chakraboty, 
Jundong Huang and Chunhui Luo for the great physics discussions we have had and also 
for their friendship. 

Finally I would like to thank my beautiful fiancee Use Friberg and my parents Avni 
and Fatmire Kryemadhi for their encouragement and support in the process of making this 
thesis happen, they have been happy for my achievements and patient with my frustrations. 



vi 



Abstract 



This dissertation describes two different projects from two different experiments. We 
have performed a search for CPT violation in neutral charm meson oscillations using data 
from the FOCUS Experiment. While flavor mixing in the charm sector is predicted to 
be small in the Standard Model, it is still possible to investigate CPT violation through 
a study of the proper time dependence of a CPT asymmetry in right-sign decay rates 
for D° — > K~ir + and D° — ► K + tt~. This asymmetry is related to the CPT violating 
complex parameter £ and the mixing parameters x and y: Aqpt °t P e £y ~~ Im£x. We 
determine a 95% confidence level limit of —0.0068 < Re£y — Im£x < 0.0234. Within the 
framework of the Standard Model Extension incorporating general CPT violation, we also 
find 95% confidence level limits for the expressions involving coefficients of Lorentz violation 
of (-2.8 < N(x,y,S)(Aa + 0.6Aa z ) < 4.8) x 10~ 16 GeV, (-7.0 < N(x, y, 5)Aa x < 
3.8) x 10~ 16 GeV, and (-7.0 < N(x, y, 6) Aa Y < 3.8) x lO" 16 GeV, where N(x,y,5) is 
a normalization factor that incorporates mixing parameters x, y and the doubly Cabibbo 
suppressed to Cabibbo favored relative strong phase S. 

We also present measurements of the A° lifetime in the exclusive decay channel A° — > J j ipA 
with J/ijj — > and A — > pit", the B° lifetime in the decay B® — > J/ipKg with J/ijj — > 

and Kg — > tt + tt~, and the ratio of these lifetimes. The analysis is based on approxi- 
mately 225 pb -1 of data recorded with the D0 detector in pp collisions at ^fs = 1.96 TeV. 
The lifetime is determined to be r(Ag) = 1.36 ± 0.30 (stat) ± 0.07(syst) ps, the B° 
lifetime r(B d ) = 1.43 ±0.12 (stat) ±0.04 (syst) ps, and the ratio 4^4 = 0.95 ± 0.22 ± 0.05. 
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In contrast with previous measurements using semileptonic decays, this is the first deter- 
mination of the A° lifetime based on a fully reconstructed decay channel. 
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Chapter 1 



The Standard Model of Particle Physics 

The Standard Model [1] is the currently accepted theory of the fundamental particles and 
their interactions. It has proven to be a successful theory. Even precision measurements 
have found no deviations so far from its predictions, with the exception of the neutrino 
masses [2]. However a number of problems remain, The Standard Model makes no room 
for including gravity and has the unattractive feature of 18 free parameters or 21 when 
incorporating massive neutrinos. And we have not measured all the parameters to the same 
precision. Also we have yet to discover the Higgs which is the cornerstone of the Standard 
Model. This short chapter will simply list its main components. 

Within the Standard Model, there are two broad categories of particles. The fundamen- 
tal fermions (fermions are particles with fractional spin) are considered to be the matter, 
the stuff of the universe. For example, the quarks in the protons and neutrons in an atomic 
nucleus and the electrons surrounding it belong to this category. The fundamental bosons 
(bosons are particles with integer spin), on the other hand, are responsible for the forces 
between the matter particles. For example, the quarks in protons and neutrons of the nu- 
cleus are held together by gluons and electrons are bound to the nucleus by the exchange 
of "virtual" photons. 

The interactions between matter particles take place through the exchange of the funda- 
mental bosons. These interactions are described by the Lagrangian of the Standard Model, 
leading to equations that specify rules for calculating quantities such as the probabilities 
for certain reactions to occur, referred to as cross sections. 
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1.1. The Fundamental Particles 
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1.1 The Fundamental Particles 

According to the Standard Model, there are 24 fundamental matter particles (see Ta- 
bles ll.ll and ll.2JI - six quarks (the up, down, strange, charm, bottom and top quarks) coming 
in three different colors [3] and six leptons (the electron, muon, tau, and a neutrino) [29]. 
All of the elements of the periodic table can be built from combinations of only three of 
these 24: the up quark, the down quark, and the electron. An oxygen atom, for example, 
has eight electrons surrounding a nucleus of eight protons and eight neutrons. Protons are 
built from two up quarks and one down quark; and neutrons are two down quarks and one 
up quark. 

Because of the strength of the forces between them (the strong force) , quarks are confined 
to exist in composites. A combination of a quark and an antiquark is a meson (e.g., a pion 
or kaon), and a combination of three quarks is a baryon (e.g., a proton or neutron). Mesons 
and baryons are collectively termed hadrons. While the vast majority of the matter we 
encounter in everyday life consists only of up and down quarks, the other four quarks are 
equally as fundamental. The essential difference is only in their greater masses. 

The lepton category of the fundamental particles consists of three negatively charged 
particles, of which the electron is prototypical, and their three very weakly interacting 
neutral partners, the neutrinos. The muon and tau differ from the electron only in mass. 
Neutrinos are only emitted during weak processes and only interact weakly, and are there- 
fore extraordinarily difficult to detect. All of these fundamental fermions have anti- matter 
partners. 

1.2 The Fundamental Forces 

The fundamental particles interact through four different forces. Probably the most 
familiar of the forces is gravity, the force of attraction between massive bodies. However 
due to the smallness of the masses of the fundamental particles, the force of gravity between 
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Table 1.1: Properties of the quarks. 



Name 


Symbol 


Charge 


Up Quark 


u 


+2/3 


Down Quark 


d 


-1/3 


Strange Quark 


s 


-1/3 


Charm Quark 


c 


+2/3 


Bottom Quark 


b 


-1/3 


Top Quark 


t 


+2/3 



Table 1.2: Properties of the leptons. 



Name 


Symbol 


Mass 


Charge 


Electron 


e 


0.511 MeV/c 2 


-1 


Muon 


I 1 


105.6 MeV/c 2 


-1 


Tau 


T 


1.777 GeV/c 


-1 


Electron 
Neutrino 




< 3eV/c 2 





Muon 
Neutrino 


v l* 


< 0.17MeV/c 2 





Tau 
Neutrino 




< 18MeV/c 2 
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Table 1.3: Properties of the gauge bosons. 



Name 


Symbol 


Mass 


Charge 


Photon 


7 








W Boson 


W ± 


80.4 GeV/c 2 


±1 


Z Boson 


z» 


91.19 GeV/c 2 





Gluon 


9 









any of them is negligible. It is hoped that gravity will someday be described by a theory 
unifying it to the other three forces, but the Standard Model does not incorporate the 
gravitational force. The forces of the Standard Model are all described by the exchange of 
force-carrying particles, the gauge bosons (see Table IT3|) . 

Electromagnetism is the force responsible for the repulsion between like charges, the 
attraction between unlike charges, the deflection of charged particles in magnetic fields, etc. 
In the Standard Model, the force is due to the exchange of "virtual" photons, which interact 
with any charged body. In addition, in Quantum Electrodynamics (QED), the Standard 
Model theory of electromagnetism, photons can convert to electrons and positrons, an 
electron can emit a photon, electrons and positrons can annihilate into photons, etc. 

The strong force acts only on quarks, and is due to the exchange of gluons. Similiar to 
the electric charge for electromagnetism, the strong force proceeds through a charge of its 
own, the "color" charge. But unlike electromagnetism, where the photon has no electric 
charge of its own, gluons do carry color charge, allowing them to interact among themselves 
and thus creating a much more complex situation. The strong force binds quarks tightly 
into hadrons, so tightly that the quarks never appear unbound. The Standard Model theory 
of the strong force is Quantum Chromodynamics (QCD). 

Finally, the weak force affects all of the fundamental particles. It is carried by the 
and Z° bosons. Nuclear beta decay is the most familiar example of this force, where one of 
the down quarks of a neutron converts to an up quark by emitting a W ~ , which subsequently 
decays to an electron and an electron antineutrino. The Standard Model theory of the weak 
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force and QED are united into a single theory, the electroweak theory [4], by introducing 
a Higgs Boson. The search for the Higgs Boson, the last of the Standard Model particles 
to be experimentally undiscovered, is one of the major efforts of contemporary high energy 
physics. 

1.3 Lagrangians 

The mathematical structure of the Standard Model is contained in a series of La- 
grangians. 

£-SM = £-QED + £-QCD + £\Veak- (1-1) 

For example, the QED Lagrangian can be written as: 

Cqed = - qAff ~ m)4> - l -F^ u F^ , (1.2) 

where 



Here, ip represents a particle with charge q, and represents the photonic vector field. 
The first term describes the kinetic energy of the particle, and together with the mass term, 
they constitue the Lagrangian density of a free particle. The second term describes the 
interaction of the particle with the electromagnetic field. The last term describes the free 
electromagnetic field. 
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Part I 



Chapter 2 



CPT Formalism 



The combined symmetry of charge conjugation (C), parity (P), and time reversal (T) 
is believed to be respected by all local, point-like, Lorentz covariant field theories, such 
as the Standard Model we outlined in Chapter ^ However, extensions to the Standard 
Model based on string theories do not necessarily require CPT invariance, and observable 
effects at low-energies may be within reach of experiments studying flavor oscillations [5]. A 
parametrization [6] in which CPT and T violating parameters appear has been developed, 
which allows experimental investigation in many physical systems including atomic systems, 
Penning traps, and neutral meson systems [7]. Using this parameterization we present the 
first experimental search for CPT violation in the charm meson system. 

2.1 Mixing Formalism 

Before we introduce CPT formalism, we want to outline mixing formalism, since much 
of the notation from standard mixing in charm mesons is used there. 

Assuming CP conservation in the charm meson system, the CP eigenstates of the neutral 
D meson can be written as, 



If we define CP\D°) = \D°), it then follows that \D{) is a CP-even state and \D 2 ) is CP-odd. 

7 




and 




(2.1) 
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The time evolution of the \D\) and I-D2) states is given by 



\Di(t)) = e-^-^^lAW), 



(2.2) 



where Mi and T{ are the mass and the width for state i. Rearranging Eq. I2.1( we find in 
terms of D\ and D 2 , that a pure D state produced at time t = is 



1 



\D°(t = 0)) = -j=(\D 1 ) + \D 2 )). 



(2.3) 



We obtain the time evolution of D° by plugging in the time evolution of the D\ and D 2 
states as given by Equation 12.21 



\D°(t)} = (V^-^'IZMO)) + e- l ^ h - iE T^\D2(0)}^ . 



(2.4) 



This can be expressed in terms of the \D°) and \D°) by using the relations in the Equa- 
tions E^and combining like terms: 



\D°(t)) = -(A + \D ) + A-\D°) 



with 



(2.5) 



A± = e -i{Mi-i^)t ± e -i(M 2 -i^)t_ 



(2.6) 



The terms A± can be arranged in more convenient forms by using the definitions: 

r _ £1 + r 2 M _ m 1 + m 2 ^_ -(M1-M2) and y _ r 1 -r 2 



2F 



The expressions for A± with these definitions are 



A, = 2e~ (iM+ 2)*cosh 



(y + ix) 



and 



A_ = -2e-( iM +?) t sinh 



(y + ix 



Tt 



(2.7) 



(2.8 



(2.9) 
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Now we pose the question: what is the probability of an originally pure D° state to decay 
to K + ir~? Define (/| as the vector representing the final state K + ir~. The amplitude for 
this decay process is: 



</l^(*)> = 
e _ (iM+ r )t f h 



(y + ix 



Tt 



(f\D°) -sinh 



(y + ix 



Tt 



</l^°>}, 



(2.10) 



where (f\D°) is the Double Cabbibo Suppressed (DCS) decay amplitude and (f\D°) is the 
Cabbibo Favored (CF) amplitude. The DCS to CF amplitude ratio is written as 

(f\D°) 



(2.11) 



(f\D Q ) 

where Rdcs is the DCS to CF branching ratio and 5 is a strong force phase between DCS 
and CF amplitudes. Plugging this in and approximating the hyperbolic functions with the 
first term of their Taylor series expansions we obtain 



(f\D°(t))=e-( m+ ^ t (f\D°) 



V RDese tS - {y + ix)^ 



Tt 



(2.12) 



Finally the probability is the absolute value square of the amplitude: 



\(f\D°(t))\ 2 = 



-n 



\{f\D°)f 



Rdcs + VRDCs(ycos5 - xsinS)Tt + ( 



(2.13) 



We use the soft pion from the decay D* + — > D°ir + to tag the flavor of the D at production, 
and the kaon charge in the decay D° — > K~ir + to tag the D flavor at the time of decay. 
Right-sign signal (RS) is obtained by requring that the soft pion charge is equal the opposite 
of kaon charge. Wrong-sign signal (WS) is obtained by requring that the soft pion charge 
is equal the kaon charge. Define a quantity R(t), which is the time-dependent rate for the 
WS process relative to CF (RS) branching fraction or 



\(f\D°)\ 2 



(2.14) 
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and define the parameters x' and y' that are related to the mixing parameters, x and y, by 
a strong phase rotation: 

x' = xcos5 + ysm5, y = ycose) — xsinJ. (2-15) 

Redefine t in units of D lifetime ( where Tt = t/r D o ) to obtain an expression for the 
lifetime evolution of the decay D° — > K + ir~: 

(2.16) 

The first term in Ea. l2.16l is a pure DCS decay amplitude, the second term is the interference 
of DCS and mixing, and the third term is a pure mixing term. 



R{t) 



Rdcs + V RDcsy't + 



,x' 2 +y' 2 



)t 2 



2.2 Proper Time Asymmetry 

The time evolution of a neutral-meson state is governed by a 2 x 2 effective Hamiltonian 
matrix A in the Schrodinger equation. For a complex 2x2 matrix, it is possible to write the 
two diagonal elements as the sum and difference of two complex numbers. It is also possible 
to write the off-diagonal elements as the product and ratio of two complex numbers. Using 
these two facts, which ultimately permit the clean representation of T and CPT- violating 
quantities, a general expression for A can be taken as [6]: 

1 fu + i vw~ l \ 

A = -AA , (2.17) 

2 \VW U-i) 

where the parameters U, V, W, and £ are complex. The requirements that the trace of the 
matrix is trA = A and that the determinant is detA = A1A2 impose the identifications 
U = A/ A A, V = yl — £ 2 on the complex parameters U and V. The free parameters in 
Eq. 12.171 are therefore W and £. These can be regarded as four independent real quantities: 
W = we luJ , and £ = Re£ + ilm£. One of these four real numbers, the argument u, is 
arbitrary and physically irrelevant. The other three are physical. The modulus of W 
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controls T violation with w = 1 if and only if T is preserved. The two remaining real 
numbers, Re£ and Im£, control CPT violation and both are zero if and only if CPT is 
preserved. The quantities w and £ can be expressed in terms of the components of A 



as w = >y | A21/A12 1 , £ = AA/AA, where AA = An — A22 and AA is the difference in the 
eigenvalues. Ai = Mi — \iT\, A2 = M2 — and AA = T(x — i y). ^ is phenomenologically 
introduced and therefore independent of the model. Indirect CPT violation occurs if and 
only if the difference of diagonal elements of A is nonzero. To determine the time-dependent 
decay amplitudes and probabilities, it is useful to obtain an explicit expression for the time 
evolution of the neutral D meson. Doing the same exercise as in the mixing formalism (with 
A as the matrix) we obtain, 



\D°(t,i,p)J \SVW- 1 C-S£j\D°J 



The functions C and S depend on the meson proper time t, sidereal time t and are given 
by: 

C = cos (\&\t) e~^ iXt S = -ism ( -AXt] e~^ xt . (2.19) 



2 



One can easily extract time-dependent decay probabibilities by manipulating Eq. 12.181 
For the decay of D° to a right-sign final state / (which could be a semileptonic mode, 
or a Cabibbo favored hadronic mode (with DCS negligible)), the time-dependent decay 
probability is 



Pf(t) EE |</|T|D°(t))| s 
■\F\ 2 eM~l 

-2 Re £ sinhA7t/2 - 2 Im £ sin Ami]. (2.20) 



i|F| 2 exp(-|i) x [(1 + |^| 2 )coshA 7 t/2 + (1 - |£| 2 )cosAmt 



The time-dependent probability for the decay of D to a right-sign final state /, Pj(t), 
may be obtained by replacing in the above equation £ — ► — £ and F — > F. In the formula, 
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F represents the basic transition amplitude for the decay D° — > /, A7 and Am are the 
differences in physical decay widths and masses for the propagating eigenstates and can be 
related to the usual mixing parameters x = AM/T and y = Ar/2l\ The complex param- 
eter £ controls the CPT violation and is seen to modify the shape of the time dependent 
decay probabilities. Expressions for wrong-sign decay probabilities involve both CPT and 
T violation parameters that scale the probabilities, leaving the shape unchanged. Using 
only right-sign decay modes, the following asymmetry can be formed, 

P 7 (t) - P(t) 

which is sensitive to the CPT violating parameter £: 

2 Re £ sinhA7t/2 + 2 Im £ sinAmi 
CPT = (1 + |£| 2 )coshA 7 £/2+ (1 - |e| 2 )cosAmt' (2 ' 22) 

We can gain insight into the anticipated experimental sensitivity by plotting these func- 
tions with some reasonable assumptions. We use 95% confidence level (C.L.) upper bounds 
on the mixing parameters x and y of 5%, which is at the upper range of the current exper- 
imental sensitivity, as discussed previously. In Fig. 12. If a) we plot the proper time decay 
probabilities for D° decay under the assumption of CPT violation at the level of Re £ = 
5% and Im £ = 5%, which are independent parameters in the framework. One sees a CPT 
violation-induced wrong-sign contribution that vanishes both at zero proper time and at 
long proper times. This causes a distortion from a purely exponential decay of a D° (and 
D°), which is then visible in the asymmetry plot, Aqpt as shown in Fig. l2.lf b). Because of 
the small oscillation frequency and short lifetime, one sees only the start of the oscillation, 
growing beyond 0.3% at long proper times. Evident from Eon. I2.2"2l is that positive values 
of Re £ and Im £ work to oppose one another in the asymmetry in a linear fashion. This is 
shown in the nearly linear behavior of Acpt m Figs. I2.lf c.d) with parameters Im £ = 5%, 
Re £ = and Im £ = 0, Re £ = 5% respectively, and consequently CPT asymmetries larger 
by a factor of 10 at long proper times. In practice, experiments will be sensitive to either 
Re £ or Im £, but not both simultaneously. 
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proper time (ps) proper time (ps) 

Figure 2.1: (a) Proper time decay probabilities for D° decay, right-sign (solid) and wrong- 
sign (dashed) and (b) A C pt with Re £ = 5%, Im £ = 5%. (c) A C p T with Re £ = 0, Im £ 
= 5%. (d) A CPT with Re £ = 5%, Im £ = 0. 

2.3 Double Cabbibo Suppressed Interference 

In the previous section, we generated the asymmetry by assuming the basic transition 
amplitudes to be (f\T\D°) = F, (f\T\D°) = 0, (f\T\D°) = F, (f\T\D°) = 0. This is fine 
as long as we are dealing with a semileptonic decay or if we neglect the Double Cabbibo 
Suppressed decays (DCS) of hadronic modes. The general framework developed in Ref. [6] 
assumes that the DCS effects are negligible. However in our forming the asymmetry we 
had to deal with DCS that were comparable to Cabbibo favored decays. We therefore will 
report the asymmetry in the previous section with DCS effects included. For the decay 
D° -► K-TT+ the basic transition amplitudes are then (f\T\D°) = F, (f\T\D°) = F DCS , 
(f\T\D°)=F, and (f\T\D°) = Fh^. 

For us to see if we can neglect DCS, we have to start with the more general framework 
where the DCS interference with the Cabbibo favored decays and the DCS term are not 
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neglected. Now we assume that (f\T\D°) = Fdcs and (f\T\D°) = -Fdcs are not neglected. 
To touch base with the formalism already developed for mixing, we will take the ratio 
(f\T\D°)/(f\T\D°) = -VR^ e- i5 and (f\T\D°)/(f\T\D°) = -^R^ e* 5 . With these 
in mind, the time dependent probability into a right-sign decay as in Eq. 12.201 including 
DCS presence takes this more general form: 



P f (t) = \(f\T\D°(t))\ 2 

= -\F\ 2 exp(-J-t) x [(1 + |£| 2 )coshA 7 i/2 + (1 - |£| 2 )cosAmt 



-2 Re £ sinhA7t/2 - 2 Im £ sin Ami - 2 y/R D cs x 

[(— sinhA7f/2 cos<5 — sinAmt sin<5) 

+(coshA7t/2 — cos Ami) x (Re£ cos5 + Im£ sin£)] 

+Rdcs x (coshA7t/2 - eosAmi)]. (2.23) 



The time-dependent probability for the decay of D° to a right-sign final state /, Dj(t), may 
be obtained by replacing £ — ► — £ and F — > F in the above equation. 

Rdcs is the ratio of DCS decay to the Cabbibo favored decay. 5 is the strong mixing 
phase between DCS decay and Cabbibo favored. Using only right-sign decay modes, we 
form the asymmetry as in Eg. 12.211 that is sensitive to the CPT- violating parameter £. In 



the case of negligible contributions from DCS decay, \(f\T\D°(t))\ = 0, we can form an 
identical asymmetry contribution as given by Eq. 12.221 



2 Re £ sinhA7i/2 + 2 Im £ sinAmt 



AcPT ^ (1 + |e|2) C oshA 7 t/2 + (1 - |e| 2 )cosAmt ' (2 " 24) 



We can form the same asymmetry as in Eq. 12.221 by using, instead of probabilities given 
by Eq. 12.201 the probabilities given by Eq. 12.231 1 . We have this additional interference term 



1 We have assumed that T is not violated and thus the phase related to T is zero, i.e., the only phase that 
enters in the probabilities due to interference is the strong phase. 
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in the Acpt expression, which ignores the small contributions in the denominator: 



A l Q PT (t) 



2 VRdcs x (coshA7t/2 - cos Ami) x (Re£ cos<5 + Im£ sine)) 
(1 + |£| 2 )coshA 7 i/2 + (1 - |£| 2 )cosAmi ' 



(2.25) 



With these new modifications, the total Acpt is the sum of the two contributions, 
Acpt = A CPT (t) + A™p T (t). From this general expression of Acpt, two different ap- 
proaches can be taken. In the first approach we assume that Eq. 15 of [6] does not hold 
but Eq. 21 of [6] is valid since it is phenomenologically introduced and is not dependent 
on the model. Equation 15 of [6] is Re £ = x Im£/y in terms of x,y mixing values. With 
small values of mixing, and the fact that D° has a relatively short lifetime, the following 
approximation is valid: 



Acpt = #e<£ y tjr - Im£ x t/r + 0(^Ir DC s){x 2 + y 2 )). 



(2.26) 



In the second approach, we assume that Eq. 15 of [6] does hold, so we have Re £ = 
x Im£/y. Under this constraint, and with small values of x,y mixing, our expression of 
Acpt takes this form 



A c 



PT 



Re£(x 2 + ?/ 2 )(t/T)- 
2x 



-^-(t/r) + V -Rdcs (x cos 5 + y sin 8) 



(2.27) 



This scenario is more strict since it requires that Eq. 15 of [6] is valid. While in the first 
approach we neglected the DCS decay interference term as small, we can not neglect it in 
the second approach because DCS interference term plays a comparable role in the Acpt- 
In this thesis and the resulting journal result [18], we consider both these scenarios since 
the question of the constraint given by Eq. 15 of [6] is an open question. 
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2.4 Lorentz Violating Parameters 

In the CPT and Lorentz- violating extension (SME) to the Standard Model [19], the CPT 
violating parameters may depend on lab momentum, orientation, and sidereal time [6,20]. 
It can be shown that [20] 

AA « ^Aa M , (2.28) 

where (5^ = 7(1,/?) is the four-velocity of the D meson in the observer frame. The ef- 
fect of Lorentz and CPT violation in the SME appears in Eq. 12.281 via the factor Aa M = 
r qi a fi ~ r <32 a /?> where a^ 1 and affi are CPT- and Lorentz- violating coupling coefficients for 
the two valence quarks in the D meson, and where r qi and r q2 are quantities resulting from 
quark-binding and normalization effects. The coefficients a^ 1 and afft for Lorentz and CPT- 
violation have mass dimension one and emerge from terms in the Lagrangian for the SME 
of the form —af l q r y^q, where q specifies the quark flavor. A significant consequence of the 
four-momentum dependence arises from the rotation of the Earth relative to the constant 
vector Aa. This leades to sidereal variations for CPT violating parameter £. In the case of 
FOCUS, a forward, fixed-target spectrometer, the £ parameter assumes the following form: 

£(t,p) = -^-[Aa + PAa z cosx + Psmx(Aa Y smQt + Aa x cosnt)j, (2.29) 
Aa 

where and t are the sidereal frequency and time respectively, and X, Y, Z are non-rotating 
coordinates with Z aligned with the Earth's rotation axis. Aax, Aay, Aclq and Aaz are 
the differences in the Lorentz-violating coupling coefficients between valence quarks, x is 
the angle between the D° momentum and the Z axis. The parameter A A in terms of x and 
V is: 

AA = -Am - ^ = xT- iyT = T(x-iy)= ^ - ^ , (2.30) 
2 r 

where r is the mean lifetime of the D° meson. 
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2.5 Sidereal Time 

Since CPT-violating parameters depend on sidereal time, we outline what sidereal time 
is and how we calculate it. Sidereal time is time according to the stars and not the sun. 
Careful observation of the sky will show that any specific star will cross directly overhead 
(on the meridian) about four minutes earlier every day. In other words, the day according 
to the stars (the sidereal day) is about four minutes shorter than the day according to the 
sun (the solar day) . If we measure a day from noon to noon - from when the sun crosses 
the meridian (directly overhead) to when the sun crosses the meridian - again we will find 
the average solar day is about 24 hours. If we measure the day according to a particular 
star from when that star crosses the meridian to when that star crosses the meridian again 
we will find the average sidereal day is 23 hours and 56 minutes long. The 4 minute lag is 
explained by the fact that the earth not only rotates about its axis but also proceeds along 
its orbit around the sun, while with respect to the stars the earth's motion around the sun 
can be neglected. Denote t as Greenwich Mean Sidereal Time (GMST), and d the number 
of days that have passed or have to be passed since the epoch. In most astronomy books, 
the epoch starts on January 1st, 2000 AD, 12:00 noon Greenwich London time. There is 
an algorithm that finds the total number of days since that epoch [21]. Let d be the total 
number of full days that have passed or are to pass since that epoch, y the number of years, 
and m the number of months: 

7 rn —I— Q 2T57TZ' 

d = 367y - int(-(y + int(— j^— ))) + int(— — ) + day - 730531.5. (2.31) 

Now we want to have the total number of full days plus the fraction of a day, so we have 
dtot = d+(h + mm/60) /24 in order to find sidereal time at that particular time of the day. 
h is the Greenwich UT hour and min is Greenwich UT minute. Now that we have the total 
number of days including the fractional part, the GMST angle [21] is given by: 



GMST ang i e = 280.46061837 + 360.98564736629 x d, 



(2.32) 
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From this number we remove multiples of 360, and what is left is the hour angle. In order 
to convert to hours we take t = GMSThour = GMST ang i e 24/360. In our experiment, spill 
number comes with a time stamp, so we had to map spills with time stamps known to within 
1 minute. A spill results in collisions and events. These time stamps were in Chicago time, 
and had to be converted to Universal Time (UT). During the year that data was taken, two 
important dates had to be included, Daylight Saving Time ending October 26, 1996 00:00 
and Daylight Saving Time begining April 6, 1997 00:00. During Daylight Saving Time, we 
have to add 5 hours to Chicago local time to find Universal Time, but when there is no 
Daylight Saving Time, we add 6 hours to Chicago Local Time. After they are converted to 
UT, Eg. 12.311 and Eg. I2.32l were used to find the Greenwich Mean Sidereal Hour t. 

2.6 Previous Searches 

Searches for CPT violation have been made in the neutral kaon system. Using an 
earlier CPT formalism [8], KTeV reported a bound on the CPT figure of merit tk = 
\mj(0 — mjco\/ m K° < (4-5 db 3) x 10~ 19 [9]. A more recent analysis, using the framework 
described in reference [6] and more data extracted limits on the coefficients for Lorentz 
violation of Aax,Aay < 9.2 x 10~ 22 GeV [10]. CPT tests in B° meson decay have been 
made by OPAL at LEP [11], and by Belle at KEK which has recently reported r B = 
\m B o — m-g \/m B o < 1.6 x 10~ 14 [12]. 

To date, no experimental search for CPT violation has been made in the charm quark 
sector. This is due in part to the expected suppression of D° — D° oscillations in the 
"Standard Model" , and the lack of a strong mixing signal in the experimental data. 

Recent mixing searches include a study of lifetime differences between charge-parity 
(CP) eigenstates from FOCUS, which reported [13-15] a value for the parameter yep = 
(3.42 ± 1.39 ± 0.74)%. The CLEO Collaboration has reported 95% confidence level bounds 
on mixing parameters x' and y' (related to the usual parameters x and y by a strong phase 
shift): [16] (l/2)x' 2 < 0.041% and -5.8% < y' < 1%. FOCUS has reported [17] a study of 
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the doubly Cabbibo suppressed ratio (Rdcs) for the decay D° — > K + ir~ and has extracted 
a contour limit on y' (of order few %) under varying assumptions of Rdcs and x'. The 
question arises - what can be learned about indirect CPT violation given the apparent 
smallness of mixing in the charm system? It turns out that even in the absence of a strong 
mixing signal one can still infer the level of CPT violation sensitivity through study of the 
time dependence of D° decays, which we show in this thesis. 



Chapter 3 



The E831/FOCUS Experiment at 
Fermilab 

FOCUS is a high-energy photoproduction experiment that took data during the Fer- 
milab 1996-1997 fixed-target run. A bremsstrahlung-generated photon beam with energies 
ranging from approximately 20 to 300 GeV was incident on a BeO target. While the pri- 
mary purpose of the FOCUS experiment is to study the photoproduction of charm and the 
properties of charmed mesons and baryons, the experiment has also been able to collect 
an impressive sample of light quark events. This chapter will give a brief overview of the 
FOCUS experiment and its detector. 



3.1 Physics Overview 



The FOCUS experiment has been at the forefront of charm physics since its analysis 
efforts began around 1998. Improving on its predecessor, E687 [22], FOCUS has been 
able to reconstruct over one million D mesons (see Fig. 13.1(1 . Over thirty papers have 
been published on topics such as semileptonic charm decays, charmed baryon lifetimes, CP 
violation in the charm sector, and the spectroscopy of charmed meson and baryon excited 
states 1 . 

A list of publications and more detail concerning ongoing physics analyses can be found at http:/ /www- 
focus. fnal.gov/ 
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Figure 3.1: The FOCUS D meson signal in three different final states. L/a > 5 was applied, 
where L/a is the Z) meson decay length over its error. 
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3.2 The Accelerator 

During the Fermilab 1996-1997 fixed-target run, 800 GeV protons from the Fermilab 
Tevatron were used to feed an array of fixed-target experiments. In the Main Switchyard, 
the proton beam extracted from the Tevatron was split into a meson beam line, a neutrino 
beam line, and a proton beam line. The FOCUS photon beam originated from the proton 
beam line. Figure l3~2*l shows the general layout of Fermilab and the fixed target lines. 

The 800 GeV protons of the Tevatron are generated in a series of five stages, each 
stage increasing the energy of the beam. The process begins with the Cockcroft- Walton, 
where electrons are added to hydrogen atoms to form negatively charged ions. The negative 
electric charge allows the H~ ions to be accelerated across an electrostatic gap to an energy 
of 750keV. Next, the H~ ions are fed into a linear accelerator (Linac). The Linac accelerates 
the H~ ions from 750 keV to 400 MeV using a series of RF cavities. Once at the end of 
the accelerator, the ions are stripped of their electrons in a thin carbon foil, the result of 
which is a 400 MeV proton beam. From the Linac, the proton beam is picked up by the 
Booster synchrotron. With a relatively small diameter of 500 feet, the Booster accelerates 
the protons from 400 MeV to 8 GeV. By way of the Main Injector, the protons are now 
ready to enter the much larger Main Ring, a synchrotron with a 4 mile circumference housed 
in the same tunnel as the Tevatron. The Main Ring brings the energy of the protons from 
8 GeV up to 150 GeV. In the final stage of acceleration, the protons are transferred from 
the Main Ring to the Tevatron. Using 1000 superconducting magnets, the Tevatron boosts 
the proton energy from 150 GeV to its final energy of 800 GeV for fixed target experiments. 

During the fixed-target run period, the Tevatron held 1000 proton bunches separated 
by 20 ns. The acceleration process went through a one minute cycle: 40 seconds were spent 
filling the Tevatron with 800 GeV protons, and then during the remaining 20 seconds the 
protons were extracted from the Tevatron and routed through the Main Switchyard. The 
FOCUS experiment was located in Wideband Hall at the end of the proton fixed target line. 
Data collection within the FOCUS experiment was divided into separate "runs," periods of 
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FERMILAB'S ACCELERATOR CHAIN 




Figure 3.2: A schematic of the layout of the Fermilab accelerator complex. 
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roughly one hour of running. 

3.3 The Photon Beam 

Once the 800 GeV protons have been extracted from the Tevatron and sent down the 
proton fixed target line, the proton beam is converted to a photon beam [23] through a series 
of stages (see Fig. 13 .3j) . The process begins 365 m upstream of the FOCUS experimental 
target where the 800 GeV proton beam interacts with the 3.6 meter long liquid deuterium 
production target. This interaction results in a spray of all varieties of charged and neutral 
particles. The charged particles are swept away by dipole magnets and collimators, leaving 
only neutral particles, primarily photons, neutrons, and K^s. These neutral particles are 
sent through a lead converter that converts most of the photons in the neutral beam to 
e + e~ pairs. The e + e~ pairs are guided around a thick beam dump using a series of dipole 
magnets, and the remaining neutral particles in the beam are absorbed by the dump. The 
series of dipole magnets leading the e + e~ pairs around the neutral particle dump consists 
of (1) the Momentum Dispersing Dipoles, the magnets that initially cause electrons to bend 
one way and positrons the other; (2) the Momentum Selecting Dipoles, which are optimized 
to select electrons and positrons with momenta around 300 GeV; and (3) the Momentum 
Recombining Dipoles, the magnets that recombine the electrons and positrons into a single 
beam. Once around the neutral beam dump, the e + e~ beam is further focused by the 
Focusing Quadrupoles. 

The e + e~ beam can now be used to generate a photon beam using the bremsstrahlung 
process. About 40 m upstream from the FOCUS experimental target, the e + e~ beam 
is sent through a lead radiator. The individual electrons and positrons radiate photons 
through bremsstrahlung. Because of the extremely high energy of the e + e~ beam of around 
300 GeV, the radiated photons travel in a direction nearly identical to the original direction 
of the e + e~ beam. After radiating, the electrons and positrons are swept into instrumented 
beam dumps (the Recoil Positron and Recoil Electron detectors) by the Sweeping Dipoles, 
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Figure 3.3: Generating the high energy photon beam used by FOCUS. 
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and only a high energy photon beam remains. The mean photon beam energy is around 
150 GeV, but in addition there is a long low energy tail reaching down to around 20 GeV. 

The energy of each photon in the beam nominally is measured by the beam tagging sys- 
tem. Before entering the Radiator, the energies of the electrons and positrons are measured 
by a set of five silicon planes interspersed between the Recombining Dipoles. After passing 
through the Radiator, when the electrons and positrons are swept to opposite sides of the 
experimental target, their energies are again measured, this time by lead-glass calorimeters, 
the Recoil Electron and Recoil Positron detectors. The energy of the radiated beam photon 
is then just the difference in energy of the electron (or positron) before and after the Radi- 
ator. In the case of a multiple bremsstrahlung event, the energy of the extra noninteracting 
photons is measured by a small central calorimeter, the Beam Gamma Monitor, and this 
energy is subtracted from the original measurement. In other words, the tagged photon 
beam energy (Ebeam) is calculated from: 

Ebeam = Ejnc — Eout — Ebgm, (3.1) 

where Ejnc is the incident electron (positron) energy before radiating, Eout is the elec- 
tron (positron) energy after radiating, and Ebgm is the energy of any additional photons 
produced in a multiple bremsstrahlung event. The energy resolution of the beam tagging 
system is approximately 16 GeV. 

3.4 The Spectrometer 

The FOCUS detector, building upon the previous E687 photoproduction experiment [22] , 
is a forward multi-particle spectrometer designed to measure the interactions of high energy 
photons on a segmented BeO target (see Fig. 13.4(1 . BeO was chosen as the target material to 
maximize the ratio of hadronic to electromagnetic interactions. The target was segmented 
into four sections to allow for a majority of charmed particles to decay outside of the target 
material. 
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Charged particles emerging from the target region are first tracked by two systems of 
silicon strip detectors. The upstream system, consisting of four planes (two stations of two 
views), is interleaved with the experimental target, while the other system lies downstream 
of the target and consists of twelve planes of microstrips arranged in three views. Once 
this initial stage of precision tracking is complete, the momentum of a charged particle 
is determined by measuring its deflections in two analysis magnets of opposite polarity 
with five stations of multiwire proportional chambers. The measured momentum is used 
in conjunction with three multicell threshold Cerenkov counters to discriminate between 
pions, kaons, and protons. 

In addition to excellent tracking and particle identification of charged particles, the 
FOCUS detector provides good reconstruction capabilities for neutral particles. Kg's are 
reconstructed using the "one-bend" approximation described in Ref. [24]. Photons and 
7T° are reconstructed using two electromagnetic calorimeters covering different regions of 
rapidity. 

Three elements of the FOCUS detector are most important for the analysis of the charm 
decay D° — ► K~tt + which we use to search for CPT violation. First, the tracking system 
provides a list of charged tracks and their momenta. Second, the particle identification 
system classifies the charged tracks as pions, kaons, or protons. Third, the triggering ele- 
ments require that events satisfy a certain number of requirements before they are recorded. 
Further information on other detector elements (e.g., the calorimeters) can be found else- 
where [25,26]. 

3.4.1 Tracking 

The purpose of the tracking system is both to reconstruct the paths particles have traveled 
through the spectrometer and to measure the momenta of these particles. The first task 
is accomplished by a series of detecting planes normal to the beam direction and placed 
at advantageous positions throughout the spectrometer. Each plane consists of an array of 
parallel silicon strips or wires, depending on the detector type, which send out a signal when 
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a charged (ionizing) particle passes through a silicon plane or close by a wire. Knowing which 
wire or strip a particle has passed near or through provides a one-dimensional coordinate 
of the position of the particle on the detecting plane. By grouping planes at various tilts, 
or views (see Fig. 13.5(1 . into stations, an (x,y) coordinate can be calculated at various 
positions of z, where z is the distance from the target, and x and y are horizontal and 
vertical coordinates, respectively. Connecting the (x, y) coordinates from station to station 
(z position to z position) results in a track, the path a charged particle has followed through 
the spectrometer. 

The second task, measuring a track's momentum, is accomplished by observing the 
deflections of the charged particle in known magnetic fields. In FOCUS, this is accomplished 
by using two different large aperture dipole magnets. The first magnet (Ml) provides a 
vertical momentum kick of 0.5 GeV/c, while the second magnet (M2) provides a larger 
vertical momentum kick of 0.85GeV/c in the opposite direction. Having different strengths 
for the two magnets allows sensitivity to a larger range of momentum. A low momentum 
track will be measured well by Ml, but may be bent out of the acceptance by M2. A high 
momentum track may not be deflected enough by Ml for a good momentum measurement, 
but will be picked up by the stronger M2. The momentum of a track is calculated by using 



as the track passes through either magnet, where Kick is the constant momentum kick of 
one of the magnets and AS is the change in vertical slope of the track as it passes through 
that magnet. 

The FOCUS tracking system consists of several distinct subsystems. The upstream 
system consists of silicon strip detectors placed among the target elements (referred to as 
the target silicon system [27]) and silicon strip detectors placed just downstream of the 
target region (referred to as the SSD system). Charged tracks are followed through the 
two dipole magnets by the downstream tracking system, which consists of five stations of 
proportional wire chambers (PWC). Three stations of PWC are between Ml and M2, and 
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Y-View X-View 




U-View V-View 

Figure 3.5: A rough sketch of the different PWC views. Charged particles passing through 
these planes ionize the gas around one of the wires, the charge is collected on the nearest 
wire, and then recorded, giving a one-dimensional coordinate of the particle. Combining 
views, the (x, y) position of a particle can be calculated at a given value of z. 
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two are downstream of M2. 

The silicon strip detectors in the upstream system are essentially reverse-biased diodes 
with charge collecting strips etched on the surface. When a charged particle passes through 
the interior of the silicon, electron-hole pairs are created. The internal electric field pulls 
the freed electrons to the surface of the silicon where they are picked up by the conducting 
strip, amplified, and registered in the data acquisition system. 

The target silicon system, the silicon strip system placed among the target elements, 
is composed of two stations of two planes of silicon strip detectors with strips oriented at 
±45° from the horizontal. The first station is between the second and third target elements, 
and the second station follows immediately after the last target element (see Fig. 13. 4|) . The 
planes are 25 x 50 mm in size (the larger dimension is vertical) , and the strips have a width 
of 25 /mi, giving 1024 different channels per plane. 

The SSD system, the second system of silicon strip detectors, begins just downstream 
of the target system and extends downstream approximately 30 cm, still upstream of the 
first dipole magnet. It consists of four stations of three planes each with the silicon strips 
oriented vertically, and ±45° from the horizontal. The stations are each 6 cm apart except 
for the last, which is separated by 12cm. The first station (i.e., the most upstream) consists 
of 25 mm long strips. In the central region the strips are 25 /mi wide and in the outer region 
the strips are 50 /mi wide. The other stations consist of 50 mm long strips, with widths of 
50 /mi in the central region and 100 /mi in the outer. 

Tracks in the upstream tracking system are found in three steps. First, clusters are 
formed within each plane. That is, regions where adjacent strips have fired are grouped 
together. By measuring the amount of charge collected, the cluster is forced to be consistent 
with having been formed by a single charged track. Second, projections are formed within 
each station. In other words, clusters within planes are joined to form a very short track 
segment within a station. Finally, tracks are formed by connecting the station projections. 
The last step is accomplished by fitting different combinations of station projections with 
straight lines and taking the best fits to be the tracks. 
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The downstream tracking system is composed of five stations of proportional wire cham- 
bers (PWC). A PWC operates on roughly the same principle as a silicon strip detector. 
When a charged particle passes through a PWC, the PWC gas is ionized and the ions drift 
through an electric field and are collected by parallel metal wires. The charge is collected at 
the end of a wire giving the one-dimensional position of a track. Arranging the PWC planes 
within a station at various tilts, or views, gives an (x, y) coordinate for a PWC station. 

Five stations of four PWC planes each are interspersed throughout the FOCUS spec- 
trometer. The planes within a station are oriented vertically, horizontally, and at ±11.3° 
from the horizontal. The first three stations (most upstream) are placed between the mag- 
nets Ml and M2, and the last two stations appear downstream of M2 on either side of the 
last Cerenkov counter (C3). The first and fourth stations have the dimensions of 76 x 127cm 
and have a wire spacing of 2 mm. The second, third, and fifth stations are 152 x 229 cm 
and have a wire spacing of 3.3 mm. 

Tracks in the downstream system are reconstructed in three steps. First, hits in the 
planes with vertical strips are connected from station to station with straight lines, referred 
to as x view tracks. The line segments formed from hits in this view are straight since it is 
the projection unaffected by the magnetic field, i.e., it is the non-bend view. Second, the 
other three views (the horizontal wires, and ±11.3° wires) are combined within each station 
to form short projections. Finally, the x-view tracks and station projections are combined 
by fitting to two straight lines, one before M2 and one after, and with a bend parameter to 
take into account the track's bending through M2. 

Once tracks have been found in the upstream and downstream tracking systems, they 
must be linked together. This is accomplished by refitting all the hits of the upstream and 
downstream tracks with three straight lines and two bend parameters corresponding to the 
amount of deflection resulting from Ml and M2. With two opportunities to measure the 
momentum, tracks can be linked by enforcing consistency. Doubly linked tracks, where 
one upstream track is linked with two downstream tracks, are allowed to accommodate the 
possibility of photons converting to e + e~ pairs that do not significantly separate until after 
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Ml. 



The momentum resolution for charged tracks depends on the momentum of the track 
and whether the track has passed through Ml and M2 or just Ml. For tracks only deflected 
by Ml, the resolution is given by: 



For low momentum tracks, the momentum resolution is limited by multiple scattering within 
the detector material. The momentum resolution for high momentum tracks is limited by 
the spacing of the wires and strips and uncertainties in the alignment of the detector planes. 

3.4.2 Particle Identification 

Particle identification in FOCUS is provided by a series of three Cerenkov counters, which 
are based on the principle that when a particle travels through a medium with a velocity 
greater than c/n, where c is the speed of light in vacuum and n is the index of refraction of 
the medium, then the particle will radiate photons. Being sensitive to these radiated pho- 
tons, a Cerenkov counter can determine whether or not the velocity of a particle is above or 
below the velocity threshold, c/n. This velocity threshold corresponds to different momenta 
thresholds for particles of different masses 2 , and this is what allows a Cerenkov counter to 
distinguish between particle types. For example, if the velocity threshold of a Cerenkov 
counter were 0.9999 c, then the momentum threshold for a pion would be 9.87GeV/c, while 
the momentum threshold for a kaon would be 34.9GeV/c. So, if a track had a momentum 

of 20 GeV/c, as determined by the tracking system, then the Cerenkov counter would fire 

2 The momentum of a particle, p, is given by p — jmv, where m is the mass of the particle, v is the 
velocity and 7 = (1 - ^)~ 1/2 . 




(3.3) 



For tracks extending through M2 the momentum resolution is: 




(3.4) 
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Table 3.1: Properties of the three FOCUS Cerenkov detectors. 



Counter 


Material 


tt^ 1 Threshold 


Threshold 


p ± Threshold 






GeV/c 


GeV/c 


GeV/c 


CI 


80% He, 20% N 2 


8.4 


29.8 


56.5 


C2 


N 2 


4.5 


16.0 


30.9 


C3 


He 


17.4 


61.8 


117.0 



if the track were a pion, but would not fire if the track were a kaon. In this particular 
example, the Cerenkov counter ideally could cleanly distinguish between pions and kaons 
for all tracks with momenta between 9.87 and 34.9 GeV/c. 

By using three different Cerenkov counters filled with gases of different indices of refrac- 
tion (see Table l3~Tjl . FOCUS can cleanly distinguish between pions, kaons, and protons over 
a wide range of momentum. Now, for example, a 20 GeV/c pion, a 20 GeV/c kaon, and a 
20GeV/c proton will all have different signatures. The pion will fire all three counters CI, 
C2, and C3; the kaon will only fire C2; and the proton will not radiate at all. Notice that 
there is ideally a clean separation between pions and kaons with momenta all the way from 
4.5GeV/c to 61.8 GeV/c. E687 used a particle identification system based only on these 
thresholds and logic tables. 

FOCUS has improved on this system by measuring the angle with which photons are 
radiated by a particle traveling with a velocity above threshold. This provides additional 
information about the particle's velocity, v, since the angle of radiation, 9, is given by 

cos6> = — . (3.5) 
nv 

Therefore the higher the velocity is above threshold, the larger the ring of the emitted 
photons. The measurement of the angle has been made possible by dividing the back of the 
Cerenkov counters into arrays of cells, with smaller cells near the center of the counter and 
larger cells further out from the center. 

FOCUS has implemented a system called CITADL for particle identification based on 
the detected rings in the counters [28]. The CITADL system works by assigning likelihoods 
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to different particle hypotheses. For example, if a particle of given momentum (measured 
by the tracking system) were a pion, then we can calculate its velocity and the angles of 
radiation and thus know which cells in which counters should have fired. The likelihood 
for the pion hypothesis is then calculated based on the status of these cells. If a given cell 
should be "on" given the pion hypothesis, and the cell was found to be "on" , then the total 
likelihood for the pion hypothesis receives a contribution of 

Lceii = (l-e^) + a-o(l-e-"), (3.6) 

where n is the expected number of photoelectrons in the cell, a is the accidental firing rate, 
and Poisson statistics has been assumed. If the cell was found to be "off", then the total 
likelihood receives a contribution of 

L cell = 1 - [(1 - e"") + a - a(l - e"")]. (3.7) 

The likelihoods are summed over all the cells in the ring of cells that should have fired given 
the pion hypothesis to give a total likelihood for the pion hypothesis: 

Ln = ^ L celh (3.8) 
cells 

Similarly, likelihoods are calculated for the e ± , K ± , and p^ 1 particle hypotheses. 

To convert the likelihoods to x 2 -like measures, the CITADL system introduces the vari- 
ables 

Wi = -2]n(Li), (3.9) 

where i indicates the hypothesis under consideration, i.e., either e + e~, tt^, K ± , or p ± . The 
Wi with the lowest value indicates the most likely particle hypothesis. Since kaons and 
pions dominate the hadronic final states, useful parameters for particle identification are 
the "pionicity", defined as 

Pionicity = W K - W v , (3.10) 
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and the "kaonicity" , denned as 

kaonicity = W W -W K . (3.11) 

Increasing the "kaonicity" requirement, for example, decreases the chances a pion will be 
misidentified kaon. 

3.4.3 Triggers 

Whenever an interesting event occurs in the detector, data must be read out and stored. 
The trigger system is responsible for discriminating between interesting and uninteresting 
events. The trigger decision takes place in several stages and there are several different 
triggers based on different physics questions. The data for the K~n + analyses included in 
the remaining chapters are obtained through the hadronic trigger. While triggering elements 
are located throughout the spectrometer and serve various purposes, the hadronic trigger 
imposes only three simple criteria on events. 

First, like all other triggers, the hadronic trigger requires a coincidence in TR1 and 
TR2. TR1 is just downstream of the target assembly, and TR2 is just downstream of the 
SSD system. A coincidence in TR1 and TR2 guarantees that at least one charged track has 
passed through the SSD system. 

Second, in addition to having tracks in the SSD system, the hadronic trigger requires 
at least two charged tracks to traverse the entire downstream tracking system. The OH 
and HxV detectors are located just after the last PWC station and are designed to count 
charged tracks. The HxV detector covers the inner region of the acceptance and the OH 
detector covers the outer region. The hadronic trigger requires either two charged tracks be 
detected by the HxV or one charged track register in the HxV and one in the OH. Both 
the HxV and the OH include a vertical gap from top to bottom to allow e + e~ pairs to 
pass. 
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Finally, a minimum hadronic energy of 18 GeV as determined by the hadronic calorime- 
ter is an additional requirement imposed by the hadronic trigger. This requirement ensures 
the presence of hadronic tracks (as opposed to e + e~ tracks). 

3.5 Data Collection 

Over the course of its running, the FOCUS experiment collected 6.5 billion events 
recorded on 5926 tapes, each tape holding 4.5 Gigabytes of data. The data was collected 
over approximately 6500 runs, each run corresponding to roughly one hour of running time. 
The data was processed in four separate stages. 

(1) PassOne was where all the major reconstruction was performed, e.g., track recon- 
struction and particle identification. 

(2) Skiml separated the PassOne output into six large superstreams based on different 
physics criteria. One of the superstreems was the hadronic meson decays also called as 
"SEZDEE" (Super EaZy DEE) stream. This stream included all hadronic decays of D 
mesons. This stream was still large (approximately 300 8mm tapes, 1.2 TB worth of data). 

(3) In Skim2, the superstreams were separated into separate substreams by requiring 
more specific physics criteria. One of the substreams of SEZDEE was tuned to select D° — > 
K~n + decays by requiring that the invariant mass of K~ir + is between 1.7 GeV and 2.1 GeV, 
have a decay length significance greater than 2.5, and a confidence level of secondary vertex 
greater than 1%. After these cuts the data was reduced to a size of « 63 GB. 

(4) In the final stage, the K~ir + data was copied to the local disks at one of the Indiana 
University High Energy Physics clusters. 



Chapter 4 



Data Analysis 

In this thesis we investigate the current experimental sensitivity for a CPT-violating 
signal using data collected by the FOCUS Collaboration during the 1996-97 fixed-target 
run at Fermilab. The analysis is also described in a journal publication [18]. 

4.1 Analysis Aproach 

The data analysis is as follows. We analyze the two right-sign hadronic decays D° — > 
K~tt + and D° — ► K + ir~ . We use the soft pion from the decay D* + — ► D°7r + to tag the 
flavor of the D at production, and the kaon charge in the decay D° — > K~7r + to tag the D 
flavor at the time of decay. (Charge conjugate modes are assumed throughout this thesis.) 

4.2 Analysis Cuts 

D° — > K _ 7T + events were selected by requiring a minimum detachment of the secondary 
(decay) vertex from the primary (production) vertex of 5(7/,. ctl is the decay length error. 
The primary vertex was found using a candidate driven vertex finder which nucleated tracks 
about a "seed" track constructed using the secondary vertex and the reconstructed D 
momentum vector. Both primary and secondary vertices were required to have confidence 
level fits of greater than 1%. The D*-tag is accomplished by requiring the D* — D° mass 
difference to be less than 3 MeV/c 2 of the nominal value [29]. 
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Kaons and pions were identified using the Cerenkov particle identification cuts. These 
cuts are based on likelihood ratios between the various stable particle hypotheses, and 
are computed for a given track from the observed firing response ("on" or "off") of all 
cells within the track's {[3 = 1) Cerenkov light cone in each of three multi-cell, threshold 
Cerenkov counters as described earlier. The product of all firing probabilities for all cells 
within the three Cerenkov cones produces a x 2 -like variable called Wj = - 2xlog(likelihood) 
where i ranges over electron, pion, kaon and proton hypotheses. For the K and the tt 
candidates, we require Wj to be no more than 4 greater than the smallest of the other three 
hypotheses (Wj — W m i n < 4) which eliminated candidates that are highly to have been 
misidentified. In addition, D° daughters must satisfy the slightly stronger Kit separation 
criteria — Wk > 0.5 for the K and Wk — W % > —2 for the tt. Doubly misidentified 
D° — > K~ir + candidates are removed by imposing a hard Cerenkov cut on the sum of the 
two separations ((W n — Wk)k + (Wk — W^)^ > 8). Primary vertices that lie in the TR1 
region are poorly reconstructed so we exclude events in TR1, by imposing the z coordinate 
of the primary vertex < 2 cm. Fig. 14.11 shows the invariant mass distribution for the two 
D*-tagged, right-sign decays D° — > K~ir + and D° — > K + ir~. Fig. 14.21 shows the invariant 
mass distributions for right-sign decays split up into particle and antiparticle. A fit to 
the mass distribution is carried out using a Gaussian function to describe the signal and a 
second-order polynomial for the background. The fit yields 17227±144 D° and 18463±151 
D° signal events. 

Fig. 14.31 shows the primary and secondary vertices for D°'s for the run period 6 which 
has 4 target segments interwoven with target silicons. Most of the primary vertices lie 
within the target segments and some in the target silicons. The contours of the target 
segments and target silicons can be seen. About 60% of D° decays occur outside of target 
segments. 

The reduced proper time is a traditional lifetime variable used in fixed-target experi- 
ments that uses the detachment between the primary and secondary vertex as the prin- 
cipal tool in reducing non-charm background. The reduced proper time is defined by 
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Figure 4.1: Invariant mass distribution for the sum of D° and TP right-sign decay candi- 
dates. 

t' = (I — Nai)/(P'jc) where t is the distance between the primary and secondary ver- 
tex, (Ji is the resolution on I, and N is the minimum detachment cut required to tag the 
charmed particle through its lifetime. Fig. 14.41 shows reduced proper time distributions for 
the two right-sign decays: D° — > K~ir + and D° — > K + ir~. 

Table l4~Tl shows a summary of the fits for D° — > K~tt + and D° — » K + ir~. It gives us 
an overall picture of yields, signal to background ratios, masses and lifetimes. 1 The above 
cuts have been chosen to maximize the signal to background ratio. 




It is useful to know how how the signal to background ratio is distributed in bins of 
reduced proper time. We denote Si as the amount of signal in bin i and Bi the amount 
of background in the same bin. When we apply sideband subtraction, each event carries a 
1 From a simple exponential fit. 
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Figure 4.2: Invariant mass of (a) D° — > i^~7r + and (b) D° — > for data (points) fitted 

with a Gaussian signal and quadratic background (solid line). The vertical dashed lines 
indicate the signal region, and the vertical dotted lines indicate the sideband region. 
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Figure 4.3: (Left) z position of the primary vertex of D's for runnumber > 9750 and (right) 
z position of -D's secondary vertices for runnumber > 9750. The dashed line is to guide the 
eye. 



Table 4.1: Summary of the fits for D° — > K ir + and D° — > i^ + 7r . 



Parameter 


D° 


D° 


D° + D° 


Yield 


18287 ± 235 


17085 ± 224 


35342 ± 322 


S/fl 


10.24 


11.11 


10.12 


Mass (MeV/c 2 ) 


1868.50 ±0.11 


1867.80 ±0.11 


1868.10 ±0.08 


d(MeV/c 2 ) 


13.51 ±0.10 


13.63 ±0.10 


13.57 ±0.07 




412.8 ±3.8 


405.7 ±3.8 


409.6 ±2.7 



weight and thus errors of each bin will depend on signal to background ratio. The smaller 
this ratio, the larger the errors. When there is only signal then the error is equal to the 
square root of the bin content. Let's see quantitatively what happens. When the sideband 
lines are chosen as in Fig. 14.21 we get a formula that connects error with signal to background 
ratio erri = y/SjJl + 1.5 x Bi/Si). Based on this formula, we extract signal to background 
ratio per each bin when we know erri and Si . Fig. 14.51 shows the distribution of signal to 
background ratio in bins of reduced proper time for D° and D°. Both show that signal 
to background ratio decreases in large t' . This is due to the fact that contamination from 
other charm mesons is more likely at larger t' values than for smaller values. 
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Figure 4.4: Background subtracted reduced proper time distributions for D° and D°. 



4.3 Results for the Asymmetry 



We plot the difference in right-sign events between D° and D° in bins of reduced proper 
time t' . The background subtracted yields of right-sign D° and D° were extracted by 
properly weighting the signal region (— 2a, +2a), the low mass sideband (—7a, —3a) and 
high mass sideband (+3<7, +7a), where a is the width of the fitted signal Gaussian. For 
each data point, these yields were used in forming the ratio: 



Y(t' 
Acprtf) = — 



'W 



Y(t') + Y(t')P } 



(4.1) 



where Y(t') and Y(t') are the yields for D° and D° and fit')-, and f(t') are their respective 
correction functions. In the absence of detector acceptance corrections, this is equivalent to 
Aqpt as defined in Ean. 12.211 The functions f(t') and f(t') account for geometrical accep- 
tance, detector and reconstruction efficiencies, and the absorption of parent and daughter 
particles in the nuclear matter of the target. The correction functions are determined using 
a detailed Monte Carlo (MC) simulation using PYTHIA [31]. The fragmentation is done 
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Figure 4.5: in bins of reduced proper time for D° (left) and D° (right). 



using the Bowler modified Lund string model. PYTHIA was tuned using many produc- 
tion parameters to match various data production variables such as charm momentum and 
primary multiplicity. The shapes of the f(t') and f(t') functions are obtained by dividing 
the reconstructed MC t 1 distribution by a pure exponential with the MC generated lifetime. 
Fig. 14.71 shows these corrections. Detector resolution effects cause less than 8% change in 
the t' distribution as measured by deviations from a pure exponential decay. The ratio of 
the correction functions, shown in Fig. I4.6f a). enters explicitly in Eg. 14. II and its effects on 
the asymmetry are less than 1.3% compared to when no corrections are applied. Due to 
the QCD production mechanism for photoproduced charm mesons, more D° than D° are 
produced in the FOCUS data sample. This has been previously investigated in photopro- 
duction by E687, in which the production asymmetries were studied in the context of a 
string fragmentation model [30]. The effect on the Aqpt distribution is to add a constant, 
production-related offset, which is accounted for in the fit. 

The ^4cpt data in Fig. 14.6( b) are fit to a line using the form of Eq. 12.261 plus a constant 
offset. The allowed fit parameters are a constant production asymmetry parameter a and 
Re£ y — Im£ x. The value of T used in the fit is taken as V = 1.6 x 10~ 12 GeV [29]. The 
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result of the fit is: 



Re£ y - Im£ x = 0.0083 ± 0.0065. 



(4.2) 



We also report a for completeness: 



a = 0.026 ± 0.009. 



(4.3) 



If one assumes mixing parameter x,y values of 5%(current 95% C.L. upper limits) and 
Im £ = 0, one obtains for Re £, Re£ = 0.17 ± 0.13. We infer one standard deviation errors 
on Re £ of 0.13, and 95% confidence level upper bounds of 0.26. 

4.4 Results for Coefficients of Lorentz Violation 

Any CPT and Lorentz violation within the Standard Model can be described by the SME 
proposed by Kostelecky et al. [19]. In quantum field theory, the CPT- violating parameter 
£ must generically depend on lab momentum, spatial orientation, and sidereal time [6,20]. 
The SME can be used to show that Lorentz violation in the D system is controlled by the 
four vector Aa^. The precession of the experiment with the earth relative to the spatial 
vector Aa would modulate the signal for CPT violation, thus making it possible to separate 
the components of Aa^. The coefficients for Lorentz violation depend on the flavor of the 
valence quark states and are model independent. In the case of FOCUS, where D° mesons 
in the lab frame are highly collimated in the forward direction and under the assumption 
that D° mesons are uncorrelated, the £ parameter assumes the following form [6] outlined 
earlier: 



iip) 



[Aao + (3Aazcosx + /3sin\;(Aaysinr2t + Aaxcosilt)]. 



(4.4) 



AA 



O and t are the sidereal frequency and time respectively, X, Y, Z are non-rotating coordinates 
with Z aligned along the Earth's rotation axis, AA = T(x — iy), and j(p) = J 1 + p 2 D o/m 2 D0 . 
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Figure 4.6: (a) The ratio of the corrections, and (b) ^4cpt as a function of reduced proper 
time. The data points represent the ^4cpt as given in Eq. 14.11 and the solid line represent 
the fit given in functional form by Eq. 12.221 (c) Re£ as a function of Greenwich Mean 
Sidereal Time (GMST). 
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Binning in sidereal time t is very useful because it provides sensitivity to components Aax 
and Aay. Since Eq. 15 of Ref. [6] translates into Re£y — Im^x = 0, setting limits on the 
coefficients of Lorentz violation requires expanding the asymmetry in Eq. 12.211 to higher 
(non-vanishing) terms. In addition, the interference term of right-sign decays with DCS 
decays must also be included since it gives a comparable contribution. One can follow the 
procedure given by equations [16] to [20] of Ref. [6] where the basic transition amplitudes 
(f\T\P°) and (f\T\P°) are not zero but are DCS amplitudes. After Taylor expansion the 
asymmetry can be written as: 



-4c 



PT 



Reg(a 2 + j/ 2 )(t/T)' 
2x 



xy 



(t/ t) + vRdcs (x cos 5 + y sin 5) 



(4.5) 



where -Rdcs is the branching ratio of DCS relative to right-sign decays and 5 is the strong 
phase between the DCS and right-sign amplitudes. We searched for a sidereal time depen- 
dence by dividing our data sample into four-hour bins in Greenwich Mean Sidereal Time 
(GMST) [21], where for each bin we repeated our fit in t' using the asymmetry given by 
Eq. 14.51 and extracted Re£. The resulting distribution, shown in Fig. I4.6f ch was fit using 
Eq. 14.41 and the results for the expressions involving coefficients of Lorentz violation in the 
SME were: 



C oz = N(x, y, 5)(Aa + 0.6 Aa z ) = (1.0 ± 1.1) x 10" lb GeV 
C x = N(x, y, 6) Aa x = (-1.6 ± 2.0) x 10" 16 GeV, 



(4.6) 



(4.7) 



and 



C Y = N(x, y, 5)Aa Y = (-1.6 ± 2.0) x 10" 16 GeV 



(4.8 



where N(x,y,5) = [xy/3 + 0.06 (x cos 5 + y sin 6)] is the normalization factor. The angle 
between the FOCUS spectrometer axis and the Earth's rotation axis is approximately % = 
53° (cosx = 0.6). We average over all D° momentum so {j(p}} ~ 7((p)) = 39 and (3 ~ 
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1. We also compare with the previous measurements for the kaon rx and B meson tb 
by constructing a similar quantity ru [8], r^> = \AA\/m D o = ^Aa^/m^o = |£||AA| = 
7(p)| Aao + 0.6 Aa^|/m D o. The result for N(x, y, S) m is: 

N(x, y, 5) r D = (2.3 ± 2.3) x 1(T 16 GeV. (4.9) 

Although it may seem natural to report rp, the parameter ru (and tk-, tb) has a serious 
defect: in quantum field theory, its value changes with the experiment. This is because it is 
a combination of the parameters Aa^ with coefficients controlled by the D° meson energy 
and direction of motion. The sensitivity would have been best if x = 90°. 



4.5 Monte Carlo 

To understand the corrections, we analyzed simulated Monte Carlo events. Our Monte 
Carlo simulation includes the PYTHIA Model for photon gluon fusion and incorporates 
a complete simulation of all detectors and trigger systems, with known multiple scattering 
and absorption effects. The default Monte Carlo flag which is responsible for scattering 
and absorption effects include a simulation of D° and D° cross sections set at half the cross 
section for a pion. The Monte Carlo was prepared such that after trigger requirement and 
analysis cuts as in the data, we reconstruct 50 times the data statistics. Fig. 14.71 shows the 
corrections fit') for D°, D°, D°+D° and the ratio of The deviations are less than 8% 
for individual f(t'). Furthermore they cancel out when we take the ratio (of the order of 
1.3%) jr&n- The ratio is the only combination used when we form the asymmetry, so our 
detector corrections on the asymmetry are very small. 

Fig. 14.81 shows the asymmetry in Monte Carlo by fitting it with the function in Eq. l2.22l 
There is enough data statistics in the Monte Carlo sample to demonstrate that there is 
no slope in the asymmetry, i.e., only a small value of Re^ = —0.0003 db 0.0210, consistent 
with zero, could result from these corrections. Thus a significant slope in observed real 
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Figure 4.7: f(t') corrections and their ratio for D° and D°. 



data should be attributed to CPT. The D°, D° production asymmetry in Monte Carlo is 
a = 0.052 ±0.001. 



4.6 Systematic Uncertainties 

Previous analyses have shown that MC absorption corrections are very small [13]. The 
interactions of pions and kaons with matter have been measured, but no equivalent data 
exists for charm particles. To check for any systematic effects associated with the fact that 
the charm particle cross section is unmeasured, we examined several variations of D° and D° 
cross sections. The standard deviation of these variations returns systematic uncertainties of 
±0.0017, ±0.3 x 10~ 16 GeV, ±0.0 x 10~ 16 GeV, and ±0.1 x 10~ 16 GeV to our measurements 
of Re£ y — Im^ x, Cqz, Cx, and Cy respectively. 

We also investigated parent (D°,D°) and daughter (K, it) absorption separately. The 
study showed that the flat corrections in MC are small, not only because absorption effects 
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Figure 4.8: Asymmetry in Monte Carlo fitted with Eq. 12.221 

are small, but also because of a cancellation due to two competing effects. The D° has a 
slightly higher absorption rate than the D°, and the net absorption rate of a (K~ , tt + ) from 
a D° is slightly lower than the net absorption rate of a (K + ,tt~) from the D°. 

In a manner similar to the S-factor method used by the Particle Data group PDG [29], we 
made eight statistically independent samples of our data to look for systematic effects. We 
split the data in four momentum ranges and two years. The split in year was done to look for 
effects associated with target geometry and reconstruction due to the addition of four silicon 
planes near the targets in January, 1997 [27]. We found no contribution to our measurements 
of Re£ y — Im£ x and Cqz- The contributions to Cx and Cy were ±1.3 x 10 -16 GeV and 
±1.6 x 10 -16 GeV respectively. We also varied the bin widths and the position of the 
sidebands to assess the validity of the background subtraction method and the stability 
of the fits. The standard deviation of these variations returns systematic uncertainties of 
±0.0012, ±0.3 x 10~ 16 GeV, ±0.9 x 10~ 16 GeV, and ±0.5 x 10" 16 GeV to our measurements 
of Re£ y — Im£ x, Cqz, Cx, and Cy respectively. Finally, to uncover any unexpected 
systematic uncertainty, we varied our l/ag and W n — Wk requirements and the standard 
deviation of these variations returns systematic uncertainties of ±0.0036, ±1.5 x 10 -16 GeV, 
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Table 4.2: Contributions to the systematic uncertainty. 



Contribut. 


Re £ y — Im £ x 


C x (GeV) 


C 0Z (GeV) 


C Y (GeV) 


Absorption 
Split sample 
Fit variant 
Cut variant 


0.0017 
0.0000 
0.0012 
0.0036 


0.0 x 10" iS 
1.3 x 10~ 16 
0.9 x 10~ 16 
1.0 x 10" 16 


0.3 x 10~ id 
0.0 x 10~ 16 
0.3 x 10~ 16 
1.5 x 10" 16 


0.1 x 10~ 1B 
1.6 x 10~ 16 
0.5 x 10~ 16 
1.1 x 10~ 16 


Total 


0.0041 


1.9 x 10 ie 


1.6 x 10 ie 


2.0 x 10 ie 



±1.0 x lO" 16 GeV, and ±1.1 x 10~ lb GeV to our measurements of Re£ y — Im£ x, Cqz, 
Cx, and Cy respectively. Contributions to the systematic uncertainty are summarized in 
Table II~2l Taking contributions to be uncorrelated, we obtain a total systematic uncertainty 
of ±0.0041 for Re£y - Im£x, ±1.6 x 10~ 16 GeV for C oz , ±1.9 x 10~ 16 GeV for C x , and 
±2.0 x 10~ 16 GeV for Cy. 



To see further details on the assessment of systematic errors, see Appendix lAl 



Chapter 5 



Conclusions 

We have performed the first search for CPT and Lorentz violation in neutral charm 
meson oscillations. We have measured: 

Re£?/-Im£2 = 0.0083 ±0.0065 (stat) ± 0.0041 (syst), (5.1) 

which leads to a 95% confidence level limit of: 

- 0.0068 < (Re f y - Im f x) < 0.0234. (5.2) 

As a specific example, assuming x = 1% and Im£ = and y = 1% (current central values 
for mixing), one finds: 

Re£ = 0.83 ± 0.65 (stat) ± 0.41 (syst), (5.3) 

with a 95% confidence level limit of 

- 0.68 < Re £ < 2.34. (5.4) 

Within the SME, we set three independent first limits on the expressions involving coeffi- 
cients of Lorentz violation of: 

(-2.8 < N(x,y,5)(Aa + 0.6Aa z ) < 4.8) x 10~ 16 GeV, (5.5) 
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(-7.0 < N(x, y, S) Aax < 3.8) x 10~ 16 GeV, (5.6) 



and 



(-7.0 < N(x, y, S) Aa Y < 3.8) x 10~ lb GeV. (5.7) 

As a specific example, assuming x = 1%, y = 1% (current central values for mixing) and 
5 = 15° (current theoretical prediction) one finds the 95% C.L. limits on the coefficients of 
Lorentz violation of: 

(-3.7 < Aa + 0.6 Aa z < 6.5) x 10" 13 GeV, (5.8) 

(-9.4 < Aa x < 5.0) x 10~ 13 GeV, (5.9) 

and 

(-9.3 < Aa Y < 5.1) x lO" 13 GeV. (5.10) 



The measured values are consistent with no significant CPT or Lorentz invariance violation. 
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Part II 



Chapter 6 



Introduction to A® 

The UA1 experiment at CERN announced the discovery of baryon in 1991 [32]. They 
measured the production fraction times branching ratio to be: 

F(A b ) x B{A b -> J/V>A) = (1.8 ± 1.0) x 10~ 3 . 

Later on, in 1996 both ALEPH [33] and DELPHI [34] measured the A£ mass in the 
decay A c 7r. Each experiment found only 4 candidates. A® was unambiguously observed 
by CDF 110 pb" 1 Run I data [35] with a mass of 5621 ± 4 (stat) ± 3 (syst) MeV, and a 
production fraction times branching ratio of 

F(A b ) x B(A 6 J/VA) = (3.7 ± 1.7 (stat) ± 0.7 (syst)) x 10" 4 . (6.1) 

Since the signal consisted of only 20 events, only a mass measurement was made; lifetime 
measurement in this mode required more data. In the second part of this thesis, we report 
a preliminary measurement of the lifetime of A® in the fully reconstructed decay mode 
A® — > J ftp A. Measuring the lifetime in this decay mode is particularly interesting, since no 
other measurement has been published in a fully reconstructed decay mode. 

More importantly, there was a long standing discrepancy in the measured value in A® and 
6-baryon lifetimes compared to theoretical predictions [29]. There was concern about this 
disagreement. Figure lo~Tl shows the lifetime ratios where the yellow bands are theoretical 
predictions. The experimental world average for r(Afe)/r(-B°) is 0.797 ±0.052, while theory 
predicted the value to be between 0.9 and 1. However the most recent calculations, only 
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in the past year, show less of a discrepancy [29]. Figure l6~2l shows the recent results. To 
understand any discrepancy, we need to consider clean decays of A® baryon such as the 
exclusive A° — ► J/ipA decay. 



T(B )/t(B u ) 



t(BJ/t(B u ) 



T(A b )/T(B°) 



x(b baryon) 

/x(B°) 



n i i i i i i i i i i i i r 



j i i I i i L 



j i i I i i L 



1.073±0.014 
1 .03 - 1 .07 

0.949±0.038 
0.99-1.01 

0.798±0.052 
0.9-1.0 

0.784±0.034 

0.9-1.0 



0.7 0.8 0.9 

Lifetime 
'orking Group 



1 1.1 1.2 
lifetime ratio 



July 2002 



Figure 6.1: World averages of b hadron lifetimes compared to older theoretical predictions 
(yellow or dark grey bands). 



The A® heavy baryon is an outstanding system to understand quark dynamics and 
Heavy Quark Effective Theory (HQET) and the operator product expansion (OPE) [39]. 
The quark content of A® is (udb) where the b quark is separated from the (ud) quarks 
that form their own spin-0 system. The A® — > J/ipA decay receives only small nonspectator 
contributions and hence its theoretical calculation is relatively straight forward. 
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Figure 6.2: World averages of b hadron lifetimes compared to more recent theoretical pre- 
dictions (yellow or dark grey bands). 



The decay mode A® — > J/ipA is said to be fully reconstructed, because all of the final 
state particles leave tracks in the detector. Thus the A® full momentum and invariant mass 
may be determined. This is in contrast to semileptonic decay modes such as Aj, — > A c iv that 
contain neutrinos. Neutrinos are neutral and interact very weakly; they leave no signal in 
the D0 detector. The semileptonic decays of A® have larger branching ratio and therefore 
are more abundant in our detector; however the sample is not as pure as in the case of fully 
reconstructed decay like A° — > J/i^A. 

Since the A^ is more massive than the B + and B° mesons, it is currently produced only 
at the Tevatron. D0 and CDF are the only currently operating detectors that can study it 
(in comparison to the B factories, Belle and BaBar, operating at the T(15)). 

The decay A° b -> J/ipA is a color-suppressed decay that proceeds through an internal 
W decay. The Feynman diagram is shown in Figure 16.31 The decay is color-suppressed 
because the colors of the quarks from the virtual W must match the colors of the c quark 
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Figure 6.3: Feynman diagram for the decay A® — ► J/i/jA. 

and the remaining diquark system. 

A brief description of the analysis follows. We reconstruct A® candidates as described in 
Chapter |H1 and first establish the A® mass signal. In addition to reconstructing A^ — > J/ipA 
we also reconstruct B% — > J/i/jK®. This serves as a control sample, since it has similar 
topology to that of A® but the B° is more abundant and has a well-known lifetime. For 
each candidate, we obtain the value of proper decay length. We then perform an unbinned 
maximum likelihood fit to the distribution of proper decay lengths in the data, to extract 
the value of ct(B°) and cr(A&). A first result was presented at Lepton-Photon conference in 
Summer 2003. The signal and lifetime presented in this thesis are found with additional data 
until January 2004. Even more data, in an analysis continued by members of CINVESTAV 
group, resulted in a submitted result [36]. 



Chapter 7 



Theory Predictions 

7.1 Spectator Model 

In the Spectator Model of hadrons, a heavy quark inside a hadron is bound to the lighter 
"spectator" quarks. When the interactions of the heavy quark with the lighter quarks are 
small, we can estimate the weak decay of the heavy quark separately. In this approximation, 
all the hadrons containing a given heavy quark have the same lifetime. This approximation 
is more valid when the quark at hand is heavier. In the simple spectator model, the decay 
width of a b hadron is: 

m _ ffiMgM ,„> 

This is the formula for the muon decay width, with the addition of CKM matrix element V c b 
for quark coupling (b to c). We assume b decays mostly to c. Plugging in nib ~ 4.2 GeV [29] 
and \V c b\ ~ 0.045 [37] this gives r = 1.2 ps. Despite the simplicity of it, the experiments 
show that this model is not sufficient. Experiments show the hierarchy of 6-quark lifetimes 
to be r(A 6 ) < t{B°) » t(B s ) < t(B+) 

7.2 Present Findings 

In the previous section we stated that the "naive" spectator model is insufficient to 
explain the hierarchy of lifetimes of heavy hadrons. In a hadron one cannot neglect the 
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strong interactions between the heavy quark and the lighter quarks, therefore a theory that 
includes these interactions is needed to explain the observed hierarchy. In the Heavy Quark 
Expansion (HQE) one uses the dimensional operators that contain l/m& terms [40]. It has 
been shown that this theory explains well the hierarchy in the £>-meson lifetimes and also 
it predicts the ratios of the lifetimes of i?-mesons. Predictions agree with experimental 
measurements. The agreement between theory and experiment gives us some confidence 
that quark-hadron duality, which states that smeared partonic amplitudes can be replaced 
by the hadronic ones, is expected to hold in inclusive decays of heavy flavors. Figure 16.21 
shows a summary of the current world average of b hadron lifetimes compared to theory. 
According to Fig. 16.21 for the B mesons, we have these experimental values and theoretical 
predictions: 

Z^A I = 1.086 ±0.017, Z^A \ h = 1.06 ±0.02, (7.2) 
r{B d y ex r{B d y th K ' 

4f4u = 0.951 ± 0.038, 414^ = L0 ° ± °- 01 ' ( 7 ' 3 ) 

which show agreement of theoretical predictions and experimental measurements. 

For a long time, the low measured value of the ratio t(A{,) /r(B d ) has been a challenge 
for the theory. According to Fig. 16.21 for the ratio r(Ai,)/T(B d ) we have: 

444U = 0.798 ±0.052, (7.4) 



r(B d ) 



\th 



0.86 ±0.05. (7.5) 



This ratio is inconsistent with the initial theoretical predictions, as shown in Fig. 16.11 
However, recent next-to-leading order (NLO) calculations of perturbative QCD [38] and 
1/irib corrections [39,40] to the spectator effects reduced the discrepancy, yielding result as 
shown in Fig. 16.21 which is in better agreement with experimental measurement. 
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7.3 Formalism 

Inclusive decay rates can be calculated in the HQE. We use the optical theorem to 
relate the decay width to the imaginary part of the matrix element of the forward scattering 
amplitude: 



F(H b ->X) = l/2m b (H b \T\H b ), T = Imi J d 4 xT{H eff (x)H eff (0)}. (7.6) 

Here H e ff represent an effective AB = 1 Hamiltonian at the scale [i = m b , 

H eff = ^V cb J2iciQ u i d ' +C2Q u 2 ' d '] + h.c, (7.7) 

where q are the Wilson coefficients, d! and u' are quark flavor eigenstates, and Q\ and Q2 
are the four-quark operators. The energy release is large in the heavy quark limit, therefore 
an OPE can be constructed for Eq. 17.61 which results in series of local operators of increasing 
dimension suppressed by powers of l/m b as shown [42]: 



r(H b - /) 



Glml , 



1927H 



\KM\ 



(gtjjiHgt) c 5 (/) (Hblbia^G^blH,) 
C3U) — — 1 o ^r-^ h 



2M Hb 



mi 



2M Hb 



+ E 



mf 2M Hh 



+ 0(l/m$) 



(7. 



where the dimensionless coefficients depend on the parton level characteristics of 

/ (such as the ratios of the final state masses to m b ). KM denotes the appropriate combi- 
nation of weak mixing angles. G^ u is the gluonic field strength tensor. The summation in 
the last term is over the four-fermion operators with different light flavors q. 

As we showed in the last section, at leading order in the heavy quark expansion, all 
heavy hadrons have the same lifetime. The situation changes at higher orders. At order 
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1/m^ the difference between meson and baryon lifetimes has to do with their structure. 
The ratio of lifetimes of A® and Bj is 

= 1 + i[/4(A&) - A{B d )] + §f[/4(A 6 ) - faBi)] + 0(l/mg), (7.9) 

where Co ~ 1.2 [42,43]. and fj,Q represent kinetic energy and chromomagnetic interaction 
corrections [42]. At this order in HQE, the difference is mainly driven by the fact that light 
quarks in A^ appear in a J p = + quantum state, reducing any correlations of spins 
between the heavy-quark and the light quark-gluon cloud. This results in ^(Aj) = 0. 
Matrix elements of kinetic energy operators cancel each other to a large degree, that results 
in a difference of at most 1-2%, which is not sufficient to explain the observed pattern of 
lifetimes. 

Dimension six operators, that enter at the 1/mf level, are the main contributors. An 
important subgroup of these operators involves four-quark operators, whose contribution is 
also enhanced due to the phase-space factor 16-7T 2 . These effects are usually called Weak 
Scattering (WS), Weak Annihilation (WA), and Pauli Interference (PI). They introduce 
major differences in the lifetimes of all heavy mesons and baryons [42-45]. Their contri- 
bution to the lifetime ratios are directed by the matrix elements of AB = four-fermion 
operators [40]: 

T S pec = T spec + T spec + T spec , (7.10) 



where Tj terms contributing to Eq. 17.61 are expressed in terms of the four-quark operators 
Of. They are defined as, 0\ = ba^(l - Ts^i^T/^ 1 ~ 7s)9j 5 0\ = 6i7 M 75^7/i(! ~ 
The recent progress has been in understanding lifetimes by concentrating on computing the 
next-to-leading order (NLO) QCD corrections to Wilson coefficients of these operators. Also 
a great deal of progress has been made in calculating matrix elements of these operators in 
quark models and on the lattice. At NLO one can parametrize the meson-baryon lifetime 
ratio as: 
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4-^ « 0.98 - (di + d 2 B)r - (d 3 ei + d 4 e 2 ) - + 4# 2 ) + <W, (7.11) 

where the scale dependent parameters di are defined in Ref. [43], and r = |^(0)| 2 /|^(0)| 2 
is the ratio of the wave functions at the origin of the A° b and B q mesons. The 8\/ m term 
represents contributions of order l/mf and higher. In Ref. [39,40], higher-order corrections 
have been calculated and they shift the ratio of the lifetimes by —4.5% in addition to 
the « 10% given by up to order l/mf. In particular, when one goes to the calculation 
of the subleading l/m& (1/mjj) as shown in Ref. [40], the inclusion of these corrections 
brings into agreement the theoretical predictions and experimental measurements of the 
ratio of lifetimes of A° baryon and Bd meson. Even though these calculations bring into 
better agreement theory with experimental measurements, one cannot conclude that the 
agreement is for certain because of the size of theoretical and experimental uncertainties. 
The discrepancy of r(Aft)/r(5d) ratio could show up again if future measurements at the 
Tevatron and later at LHC would find the mean value to stay the same, with shrinking 
errors. In this thesis, we add another clean contribution to the world-average measurement 
of the lifetime ratio. 
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The D0 Detector for Run II 

This chapter describes the D0 detector and the Tevatron upgrade at Run II. It is based 
on Diehl's review [46]. Fig. 18.11 shows an elevation view of the D0 detector. 

First in section 18.11 we define the coordinate system used in D0 detector. Section 18.21 
describes the upgrades to the world's highest energy accelerator, the Tevatron, and its 
current status. Since the tracking system and muon spectrometer are crucial to B physics, 
emphasis will be placed on the muon detector and central tracking system in this chapter. 
Section 18.31 describes the new silicon vertex system. Section 18.41 describes the upgrades to 
the central fiber tracking system. Section "8 . 51 describes upgrades to the calorimeter systems. 
Section l8~B1 describes the upgrades to the muon detectors. Section 15*771 describes the trigger 
systems. Section "8.81 describes the data acquisition system. 

8.1 The D0 Coordinate System 

At D0, the primary coordinate system has the z-axis along the proton beam direction, 
and the positive y-axis pointing up, so that (x, y, z) make a right-handed Cartesian system. 

Sometimes cylindrical (r, (f), z) coordinates are used, as are spherical (r, cj>, 9) coordinates, 
r and 8, giving respectively the perpendicular distance and the angle from the z-axis. The 
azimuthal coordinate, 0, gives the angle from the x-axis in the x — y projection. 

The angular variables are defined so that 6 = is along the positive z-axis direction, 
and = 7r/2 is parallel to the positive y-axis. 
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Figure 8.1: Elevation view of the upgraded D0 detector. 
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The rapidity y is defined as 




(8.1) 



where E is the energy, and p z is the particle momentum on the z-direction. 



The variable of pseudorapidity is often more convenient, which is defined as 



77 = — ln[tan — ]. 



(8.2) 



In the limit that m <C E (where m is the invariant mass), the pseudorapidity approxi- 
mates the true rapidity. 

"Transverse" momentum (pr) is also commonly used, pt is the momentum vector 
projected onto a plane perpendicular to the beam axis: 



This is particularly useful due to the fact that in a high energy collision, many of the 
products of the collision escape down the beam pipe, so the momenta along the beam of the 
colliding partons are unknown. However, their transverse momenta are very small compared 
to their momenta along the beam, so momentum can be considered to be conserved in the 
transverse plane. 

8.2 The Run II Accelerator Upgrade 

In Run II, beginning in March 2001, the Tevatron collides protons with antiprotons at 
a center-of-mass energy of 1.96 TeV, which is a slightly higher energy than the 1.8 TeV 
available in Run I. In the first 3 years, four times the Run I integrated luminosity has been 
collected (by the end of August 2004) . We show in Fig. 13.21 a schematic of the Fermilab 
Tevatron Collider. 



Pt = P sin 
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In Run II, a large advance has been made with the construction of the Main Injector, 
a 150 GeV synchrotron, built in a separate tunnel from the Tevatron. The Main Injector 
has replaced the function of the Main Ring for antiproton production. It produces 2 x 10 11 
antiprotons per hour [47], four times the rate of the old Main Ring. A new permanent 
magnet "Recycler Ring" [48] allows recovery and reuse of uncollided antiprotons when 
the Main Injector has produced enough to merit injecting a new store into the Tevatron. 
Instantaneous luminosities in the range 5 — 20 x 10 31 cm _2 s _1 were available early in Run II. 
Because the number of protons per bunch is near the limit of the Tevatron, the number 
of bunches is increased from Run I's 6 bunches of protons and 6 of antiprotons to 36 of 
each species. With an instantaneous luminosity of 20 x 10 31 cm _2 s _1 , an average of 5.2 
pp collisions occurs each bunch crossing. Increasing the number of bunches decreases the 
average number of collisions per bunch crossing. By decreasing the bunch width as the store 
is depleted by collisions, the luminosity can be maintained at an optimal level. 

Such changes have profound implications for the detectors. Where in the past, there 
were 3.5 microseconds between each beam crossing, in Run II this is 396 ns for 36 bunches. 
This required all the front-end electronics from the detector to be replaced with electronics 
capable of faster response. Aside from the issues of shorter times between collisions, the 
detector improves its ability to identify leptons, energetic photons and charged particles, 
and particles emerging from a secondary vertex (daughters of long-lived parents) . Increasing 
the collision rate by an order of magnitude also requires substantial upgrades of the trigger 
and data acquisition systems to handle the increased data flow. 

8.3 Silicon Microvertex Tracker (SMT) 

The purposes of the silicon detectors are to identify tracks of charged particles that decay 
before they reach the outer tracker, to reconstruct the decay vertices of b hadrons, and to 
extend or improve the track-finding and momentum resolution. Because the proton and 
antiproton beams are very thin (<7 tra nsversc < 40 fim) and because their transverse positions 
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can be maintained very close to the center of the beampipe, a typical 6-hadron decay vertex 
can be identified by the impact parameters of the tracks of ~ 300 /im. This requires fine- 
pitched detectors oriented perpendicularly to the trajectories of the charged particles and 
mounted as close as practical to the decay vertices. The fine pitch and small radius necessary 
to achieve sufficient impact parameter resolution also improves the momentum resolution 
of the tracker systems. 

The accelerator parameters contribute to the design decisions. The interactions are 
spread out over a length a z = 25 cm, setting the length scale and motivating the detectors 
to have extended geometries. The expected integrated luminosity implies a radiation dose 
of < 1 Mrad during the life of the silicon detector. This, in turn, forces the sensors and 
readout chips to be radiation-hard and operated at cold temperatures, which mitigates some 
of the adverse effects of radiation. 

The basic detector design is the "barrel geometry" comprised of layers of silicon detectors 
arranged in plates perpendicular to the beampipe. In addition to measuring separated 
vertices and improving the momentum resolution, D0 chose to use silicon detectors to 
extend the tracking coverage to high pseudorapidity with disk-shaped assemblies comprised 
of wedge-shaped silicon wafers arranged in plates oriented perpendicular to the beamline. 

The D0 silicon microvertex detector [51] consists of six barrel modules, twelve small 
disks ("F-disks"), and four large disks ("H-disks"). The mixed barrel/disk geometry pro- 
vides silicon sensors arrayed at normal incidence, as is optimal for good tracking resolution, 
to charged particles with \rj\ < 3. Figure I&21 shows the D0 silicon barrel modules with 
barrel detectors parallel to the beamline and disks perpendicular to the beamline. 

The six barrel modules are constructed with 4 layers of ladder assemblies with con- 
siderable overlap. Each barrel module is 12.4 cm in length and the total length of the 
barrel is 76.2 cm. The 4-layer coverage corresponds to the region r] = 1.5 for interactions 
at z = 0. Barrel layers one and three are constructed from double-sided silicon sensors 
with axial (r — z) and 90° stereo layers, except for the modules on each end of the barrel, 
modules one and six, which have single-sided axial strips. These are single-sided because 
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Figure 8.2: D0 silicon detector. The figure shows the configuration of the 6 barrel modules, 
the 12 "F-disks", and the 4 "H-disks". 

stereo tracking is dominated by the information from the F-disks. Barrel layers two and 
four are constructed from double-sided detectors with axial and 2° stereo layers in all six 
barrel modules. All of the detectors are AC-coupled to the readout electronics. 

The 12 "F-disks" are comprised of 12 trapezoidal wedges arranged into a plate with 
a hole for the beampipe. The active area inner radius is at 2.5 cm from the center of 
the beampipe and the outer is at 9.8 cm. The detectors are AC-coupled, double-sided, 
with strips angled at ±15° with respect to the vertical. The pitch of the p-side and n-side 
detectors is 50 /im and 62.5 //m, respectively, and the silicon wafers are 300 /xm thick. 
Naturally, the strips are of different length depending on their locations on the wedge. The 
six innermost F-disks attached to the outer sides of the six barrel modules. The three 
additional F-disks are attached to each outer side of the barrel assembly, with the effect 
of extending the acceptance of the silicon system to higher pseudorapidity, especially for 
interactions which occur at larger z. 

The four large "H-disks" are located at z = ±94 cm and z = ±126 cm. The inner radius 
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of the active area is at 9.6 cm from the center of the beampipe and the outer is at 23.6 
cm. The detectors are AC-coupled, single-sided, with 40 /xm pitch strips (pairwise readout 
makes the effective pitch 80 /im) angled at ±7.5° from the vertical. Each plane has wedges 
glued together back-to-back to provide a 15° stereo angle. These forward disks are necessary 
to provide track stubs for forward particles which would otherwise exit the region of full 
solenoidal magnetic field without hitting the outer tracker. They cover a pseudorapidity 
range of approximately 2 < \rj\ < 3. 

All of the three detector types are able to withstand a radiation dose greater than 1 
Mrad. The dose depends on the detector location and integrated luminosity. The innermost 
layer is expected to receive 0.5 Mrad/fb _1 . The effect of the radiation damage is to increase 
the leakage current, increase the bias voltage necessary for full depletion, and decrease 
the signal-to-noise ratio. The effects are temperature dependent and can be reduced by 
operated the detectors at low temperatures. The detectors are operated at temperatures 
between 5 — 10° C using a cooling mixture of deionized ethylene glycol and water. 

The D0 silicon detectors are read out using the SVX II chip [52]. SVX II is a 128 
channel, radiation hard CMOS chip mounted directly on the "High Density Interconnects" 
(HDI's), kapton-based flexible circuit, wire-bonded to the sensors at the ends of the ladders 
and wedges. Each of the 128 channels in an SVX II chip features a preamplifier, a 32 stage 
analog pipeline with 4 fj,s delay, an 8-bit Wilkinson type ADC, and a latch-based sparsified 
readout. The chip is programmable for any interaction time from 132 to 396 ns. The 
pipeline depth, ADC ramp rate, preamplifier bandwidth, and thresholds are downloadable 
to each SVX II chip. 

Connections to the outside world continue from the HDI's to 8-ft long "low-mass" cables 
that join the ends of the HDI's to unpowered "transition cards" mounted on the ends of the 
central calorimeter cryostat. Ultimately, these signals are gathered by a sequencer board 
connected by optical link to VME readout buffer electronics in the movable counting house. 
The D0 silicon detector has 792,576 channels. 
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8.4 Outer Tracker: Central Scintillating Fiber Tracker (CFT) 



The outer tracker [46] is the charged particle tracker at largest radius within the calorime- 
ter. Outer trackers perform two functions. The first is measuring the momentum and charge 
of particles produced in the collision, and the second is to provide pattern recognition as- 
sistance for the silicon detectors. 

A particle with non-zero charge q and momentum p in a solenoidal magnetic field along 
the z-direction of strength B will travel in a helix with radius r given by 



where px = y + Py- Therefore, by measuring the track's curvature in the r — <fi plane, 
we effectively measure pj>. By measuring the track's direction in the r — z plane, we measure 
Pt/Pz, which completes our measurement of the 3D momentum vector of the particle. 

Tracks in these detectors typically have several dozen hits, which allows for highly 
efficient and pure identification of these tracks. Silicon trackers in these detectors have of 
order a half-dozen hits, and while they can be used standalone, a much better way to use 
them is to find the track in the outer tracker and project this track back into the silicon. 
Once the approximate trajectory in the silicon detector has been established, the tracking 
algorithm can search the silicon detector for hits and use these hits to improve the track 
measurement. This technique uses each component to its best advantage: the silicon tracker 
measures the track's point of origin and initial direction, and the outer tracker measures 
the track's momentum. 

The D0 Solenoid and Scintillating Fiber Tracker 

One substantial improvement to the D0 detector for Run II was the addition of a supercon- 
ducting solenoidal magnet providing a solenoidal field of 2T parallel to the z-axis. Having 
a magnetic field enhances the D0 detector. It has the ability to measure the momentum 



Pt 



(8.4) 





8.4. Outer Tracker: Central Scintillating Fiber Tracker (CFT) 



72 




Silicon Microvertex Tracker 



Figure 8.3: D0's scintillating fiber tracker. 

and charge of leptons and hadrons. The muon momentum can be measured in the central 
magnetic field, before the muon has scattered in the steel of the muon toroids, and this 
allows for the comparison of the muon's momentum as measured in the solenoid to the 
momentum as measured in the toroids. We can see that the muon momentum resolution is 
much improved after matching a muon local track with central track. 

Mechanically, the D0 outer tracker is a simple device. Layers of scintillating fiber are 
placed on carbon-fiber composite cylinders. When a particle travels through one of these 
fibers, the scintillator emits light, which is totally internally reflected down the fiber, coupled 
to a clear fiber, and transported to a solid-state light detector. Fig. 18.31 shows a side view 
of the tracker. 

In detail, the D0 central fiber tracker (CFT) [53] consists of scintillating fibers mounted 
on eight concentric cylinders made from a composite of layers of high modulus of elasticity 
carbon fiber sandwiching structural foam. The cylinders range in radius from 20 to 50 cm 
and are about 2.5 m in length, except for the inner two, which are about 1.7 m in length. 
The scintillating fibers are double clad, 835 fim in diameter, and are constructed in ribbons 
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each 128 fibers wide composed of a "doublet" layer of fibers with the centers of one of 
the single layers in the space between the fibers of the other single layer. There are eight 
doublet axial (aligned along the beam axis) layers of scintillating fiber, as well as eight 
doublet stereo layers that make a ±3° angle with the beam axis. The outer (8th) layer is at 
the largest possible radius. The 7th layer is as close to the outer layer as is possible. The 
inner layer is at the least possible radius. The detector is divided into 80 sectors in cf>. Each 
pie-shaped slice has 960 fibers and the entire detector therefore has 76,800 channels. 

The scintillation light from the fibers is totally internally reflected down the length of 
the fiber. A connection is made at the edge of the tracker to a clear fiber that transports 
the light to a solid-state light detector called a Visible Light Photon Counter (VLPC). 
The number of photons available at the VLPC's, for a charged particle which traverses 
through the center of a fiber, varies from 10 to 40, depending on the charged particle's 
pseudorapidity and on attenuation due to the distance from the clear fiber connector. The 
VLPC's are small silicon devices which have an array of eight photo sensitive areas, each 
1 mm in diameter, on their surface. They operate at temperatures from 6 to 15 K, have a 
quantum efficiency of over 80%, and have a gain of 20,000 to 50,000. The high quantum 
efficiency is important because of the low number of photons. 

The momentum resolution is currently ct(pt)/pt = 0.13% • pr when scintillating fiber 
tracker measurements are combined with the silicon tracker. The pr-dependent impact 
parameter resolution is currently 30 /im in the x and y direction, and 42 fim in the z 
direction for high momentum tracks, if SMT and CFT information are combined. 

8.5 Calorimeters and Preshower Detectors 

Because of the dependence of the Tevatron's physics program on lepton identification 
and jet energy measurement, calorimetry is a critical aspect of Run II experiments. The 
calorimeters are used for identification and measurement of the electron, neutrino and jets 
from the decays of the top quarks in these events. 
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The calorimetry is divided into two parts, electromagnetic and hadronic. Nearest the 
vertex is the electromagnetic calorimeter that measures the energy of electrons and photons 
(including those from ir° and rj decay) and has improved position measurement at the 
point of maximum shower development. Farther in radius is the hadronic calorimeter which 
measures the energy of hadrons as they interact with the material of the calorimeter. Muons 
deposit a small amount of energy (due to ionization) in both sections, and the lack of a large 
energy deposit can be used to identify a particle as a muon. Neutrinos deposit no energy 
at all in these calorimeters, but the absence of energy deposition appears as a momentum 
imbalance in the transverse plane, also called "missing Ej* . 

D0's hermetic, radiation-hard uranium and liquid-argon calorimeter [50] consists of 
three separate cryostats: the Central Calorimeter (CC), and the two Endcap Calorimeters 
(EC's). Each is segmented into an electromagnetic section, a hadronic section, and a coarse 
hadronic section (inside to outside), with many layers of sampling. Each is divided into 
pseudoprojective towers covering 7?x</> = 0.1x0.1 rad. The readout of the electromagnetic 
section has four layers of longitudinal segmentation. The third electromagnetic layer, at 
EM shower maximum, has segmentation 7] x (j) = 0.05 x 0.05 rad. The readout of the 
hadronic sections have 4 (5) longitudinal layers in the CC (EC's). There are no projective 
cracks. The calorimeter provides hermetic coverage to \rj\ < 4. The energy resolution is 
ge/E = 15%/^E(GeV) © 0.4% for electrons and photons. For charged pions and jets the 
resolutions are approximately 50%/^E(GeV) and 80%/ y/E(GeV), respectively. 

The Inter Cryostat Detectors (ICD's) augment the D0 liquid argon calorimeters by 
providing a measurement of the energy in between the central and endcap cryostats. This 
improves the energy measurement for jets that straddle the intercryostat region and im- 
proves the resolution of the missing transverse energy measurement. Sixteen new ICD 
detector segments form an annular ring of 1/2" thick scintillator covering 1.1 < < 1.4 
on the hadronic section of the inner end of each EC cryostat. The sixteen segments are 
further segmented into sections of size A 77 x A0 = 0.1 x 0.1. Each section has an embedded 
wavelength shifting fiber to collect light. These are in turn connected to long clear fibers 
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which carry the light to photomultiplier tubes located underneath the cryostats, in a region 
with reduced magnetic field from the solenoid. 

The primary purpose of the Central and Forward Preshower Detectors (CPS and FPS) 
is to exploit the difference between energy loss mechanisms of electrons and photons with 
that of the backgrounds, principally hadronic jets with leading 7r°'s, to improve the trigger 
and offline purity. Secondly, they provide a precision measurement of the starting point 
of the electromagnetic showers. The CPS is cylindrically shaped, mounted on the outside 
of the solenoid magnet, and covers the region \rj\ < 1.3. The FPS are shaped like annular 
rings, mounted on the inside of the EC's, and cover the region 1.5 < \rj\ < 2.5. The FPS 
and CPS are shown in Fig. 18.41 

The D0 Run II luminosity monitor (Level-0) consists of two arrays of plastic scintillation 
counters located on the inside faces of the EC's and arranged symmetrically around the 
beampipe. The pseudorapidity coverage is 2.7 < \rj\ < 4.4. Because the solenoid field is 
~ 1 Tesla in that region, short magnetic field resistant photomultiplier tubes are used to 
read the light. A coincidence of hits in the counters on both sides of the interaction point 
provides the simplest indication that an inelastic collision occurred. 

A new near beam detector called the "Forward Proton Detector (FPD)" is available for 
Run II. It is a series of small, retractable scintillating fiber detectors placed a few millimeters 
from the beamline in the region 20-60 meters from both sides of the interaction point. They 
are triggered by small scintillation counters and read out by multi-anode PMT's. Their 
purpose is to identify scattered protons and anti-protons in diffractive events. 

8.6 Muon Systems 

The muon detection strategy at D0 relies on the penetration power of muons. Several 
meters of absorber (including the calorimeters) absorbs the vast majority of hadrons, and 
any charged particle that penetrates this material is inferred to be a muon. Because they 
are at large radius, muon detectors are large, and to keep costs reasonable, they have very 
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Figure 8.4: One quarter r — z view of the end of the D0 trackers and the start of the EC, 
indicating the Central Preshower, the Level-0 detector, the solenoid magnet and calorimeter 
cryostats, and the Forward Preshower detector. The Forward Preshower detector is shown 
in detail in the inset. 



8.6. Muon Systems 



77 



coarse granularity: typically they are single wire proportional chambers with drift times in 
excess of a microsecond. Fig. 18.51 and 18.61 show an r — z view of the D0 detector, and all 
components are illustrated in the graph. 

The muon detector consist of scintillator and drift tubes, with effectively complete cov- 
erage out to 1 7/| < 2. As seen in the layout, the detector is split at \q\ of 1 into a central and 
forward system. Each has 3 layers (usually called A,B,C with A between the calorimeter 
and iron and the other two outside the iron) of drift tubes. In the central region are pro- 
portional drift chambers (called PDTs). In the forward region are minidrift tubes (called 
MDTs). Figure IO shows the layout of PDTs and MDTs. 

In Run II, scintillator counters are added adjacent to the chambers. Scintillator has a 
response time measured in nanoseconds, so that the coincidence between the counter and 
the chamber unambiguously determines the bunch crossing. Because the timing resolution 
is substantially better than the minimum required to resolve bunch crossings, we gain 
the ability to use timing to reject certain backgrounds: particularly out of time particles 
produced upstream of the interaction point and cosmic rays. 

There is 2 or 3 layers of scintillator coverage with the forward scintillators sometimes 
called pixels, the central A-layer counters called A-phi, and the BC counters called the 
cosmic cap. Scintillator time is read out with both a 15-20 ns "trigger" gate and a 80- 
100 ns "readout" gate. The D0 detector relies on layers of scintillation counters and drift 
chambers to identify muons and measure their position and momentum. 

The muon system has three large toroid magnets, one central and two forward, which act 
as absorber in addition to the calorimeter, provide a structure on which to mount the muon 
detectors, and provide a bend to the muons so the detectors can measure the momentum. 
The calorimeters contain between 7 and 10 interaction lengths of material, depending on the 
pseudorapidity. The thickness of the calorimeters plus the toroid magnets varies between 13 
to 14 interaction lengths for |r/| < 0.9, 10 to 15 interaction lengths for 0.9 < \r]\ < 1.2, and 
18 to 20 interaction lengths for 1.2 < \rj\ < 2.0. Muons with momentum greater than ~ 1.4 
GeV/c (~ 3.5 GeV/c) penetrate the calorimeter (toroid magnet) at rj = 0. In addition, new 
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Figure 8.5: Half r — z view of the D0 muon subdetector. 
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Figure 8.6: Illustration of the D0 muon subdetector components. 




Figure 8.7: Layout of PDTs and MDTs. 
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massive shielding structures isolate the muon detectors from backgrounds generated near 
the beampipe and accelerator elements. 

8.6.1 Central Muon System Drift Chambers 

The central muon drift chambers were retained from Run I, but their electronics have been 
replaced. The drift chambers are made from extruded aluminum cells of 4-inch width and 
lengths up to 228 inches. The wires in the cells are parallel to the field in the toroid magnets 
so that the bending of the track in the toroids takes place in the drift ordinate. Refer to 

Fig, ma esi 

Individual drift chambers (PDT's) in the C- and B-layers consist of three staggered 
decks of up to 24 cells each. Drift chambers in the A-layer consist of four staggered decks of 
24 cells each, except for the ones on the bottom, which are three deck PDT's. The top and 
bottom of each cell has a copper-clad cathode pad. The copper has a milled cut separating 
it into an inner and outer pad such that the width of the inner pad alternately increases 
and decreases along its length. The wavelength of the vernier is 24 inches. Pairs of wires 
are connected through a delay chip at the end away from the front-end electronics. Fig. 18.81 
shows the geometry of a PDT cell. 

On passing through a cell, a muon will cause a hit in the cell and a hit in the neighbor 
cell which is some time later depending on the muon's proximity to the far end. Charge is 
accumulated on the inner and outer pads of the cell through which the muon passed. The 
drift time is derived from the sum of the two cells times. The distance along the wire is 
derived from the difference. The normalized difference of the integrated pad charge provides 
the distance along a pad wavelength. 

The Run I A-layer PDT's were rebuilt so as to increase their effective lifetimes. The 
ends of the PDT's were removed. The cathode pads, which outgassed a dielectric that 
coated the anode wires in Run I, were replaced with new G-10 pads that do not outgas. 
The lifetime of the B- and C-layer PDT's is long enough so that aging won't pose a problem 
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Figure 8.8: Geometry of a PDT cell. 



8.6. Muon Systems 



83 



for these detectors. 

The PDT's use a drift gas composed of 80% argon, 10% CF 4 , and 10% CH 4 . The 
maximum drift time is approximately 500 ns, longer than the bunch spacing. This poses no 
problem for the electronics, which records each hit in as many crossings in that it could have 
occurred. The time measurements are made in 1.8 ns bins. The drift ordinate resolution is 
~ 500 /im per hit, limited by the fluctuations in the drift time due to the gas. Normally 
the drift velocity is about 0.1 mm/ns. 

The muon momentum is calculated from the bend in the toroid magnet as determined 
from the difference in slopes between the line formed from the interaction point and the 
A-layer hits and the line through the B- and C-layer hits. The momentum resolution is 
expected to be a(l/p) = 0.18/p <g> 0.005 with p in GeV/c. This momentum resolution is 
worse than that expected from the fiber tracker measurement, especially at low px- 

8.6.2 Central Muon System Scintillation Counters 

An important part of the upgrade is two new layers of scintillation counters. These detectors 
not only tag the bunch crossing from which the muons originate for the slow drift chambers, 
but also reject background particles which leave hits at times other than expected from a 
muon originating at the interaction point. Refer to Fig. 18.51 18.61 

The A-layer contains 630 "A-phi" counters, each approximately 32" long and spanning 
4.5° in azimuthal angle. There are 9 counters in a row spanning —1.0 < r/ < 1.0. Each row 
has a slight overlap with a neighboring row so as to minimize the cracks between counters. 
The counters are made from 1/2 inch thick Bicron scintillator plates with many Bicron 
BCF'92 wavelength shifting fibers embedded in deep grooves. The fibers collect and trans- 
mit the light to a single photomultiplier tube (PMT). The counters have a time resolution 
of a = 4 ns and are expected to discriminate between muons produced in the collisions and 
the background, that is composed of particles backscattered from the calorimeter exit, and 
arrives 14 ns later than a muon. The "A-phi" counters span 93% of the azimuthal angle. 
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There is a gap in the A-layer coverage where the calorimeter is supported by the detector 
platform. 

A layer of scintillation counters has been added on the outside of the muon toroid 
magnet. The 240 "Cosmic Cap" C-layer scintillation counters were deployed late in the 
second half of Run I and previously have been described in detail [55]. These counters are 
between 81.5 and 113 inches long and 25 inches wide. Eight of them are mounted on the 
outside of each C-layer PDT on the top and sides of the central muon detector. Underneath 
the toroid magnet, the three layer coverage is broken up because of the support structure for 
the central platform and toroid magnets. 120 new "Cosmic Bottom" counters are arrayed 
on bottom C-layer and B-layer PDT's. The Cosmic Cap and Cosmic Bottom counters 
are made from 1/2 inch thick Bicron scintillator with Bicron BCF'91A wavelength shifter 
fibers embedded in grooves. Each counter is read out with two photomultiplier tubes. One 
advantage of using two PMT's is that coincidental tube noise is improbable. Another is the 
immediate redundancy available in case one of the PMT's fails. 

The scintillation counters have an LED pulser calibration system capable of providing 
a clocked, timed, amplitude-controlled photon pulse. Each PMT is connected by a light- 
shielded optical fiber to a light-tight box which houses bundles of LED's glued into a clear 
acrylic block. A single box may provide the photon pulse for up to 100 PMT's. The stability 
of the photon pulse is monitored with a light-sensitive diode housed in the clear block. This 
system allows the timing and amplitude to be monitored and controlled. 

The muon system has a wide range of options available for triggering. Triggers may 
be composed of coincidences of in-time hits in scintillation counters, hits in layers of the 
PDT's, and tracks found in the central tracker. Low-pr muons would rely on hits only in 
the A-layer detectors and fiber tracker. High-pr muons would also use scintillation counters 
and PDT's in the B- and C-layers. 
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Figure 8.9: D0 forward muon mini-drift tube (MDT) plane. The octant boundaries are 
shown. 

8.6.3 Forward Muon System Drift Tubes 

The Mini-Drift Tube (MDT) system [56] is comprised of three planes of drift-tubes, with 
one plane in front of, and two planes behind the forward toroid magnet. The layers are 
divided into octants with tubes of different length depending on position in the octant. As 
in the central region, the MDT A-layer has four decks of drift tubes and the B and C-layers 
have three decks each. The drift tubes, made from long aluminum extrusions of eight 1 
cm square cells, are contained in plastic sleeves. Wires in the cells are oriented parallel to 
the magnetic field of the forward toroid magnet. The sleeves of tubes are mounted on an 
aluminum support structure which also provides mechanical support for the infrastructure. 
A plane of MDT counters is shown in Fig. 18.91 

The MDT's use a non-flammable gas mixture composed of 90% CF4 and 10% CH4, with 
a maximum drift time of about 60 ns. The momentum resolution is limited by multiple 
Coulomb scattering in the iron toroid and the hit resolution of the detector. The MDT 
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Figure 8.10: Two D0 forward muon pixel octants. 

electronics uses a coarse digitization of the drift time (18.8 ns time bin). The momentum 
resolution is roughly (j{pt)/pt = 0.2. Importantly, it is on par with the resolution of D0's 
central tracker in the forward region where the full coverage of the fiber tracker's layers has 
ended. 

8.6.4 Forward Muon System Scintillation Counters 

Three layers of "Pixel" scintillation counters [56] are added to the forward region (approx- 
imately 1.0 < \r]\ < 2.0). Their primary role is muon triggering. The ~ 4800 pixel counters 
have segmentation 0.1 x 4.5° in At? x A(f>. Most of the trapezoidal shaped pixel counters 
are made from 0.5" thick Bicron scintillator with wavelength shifting bars. They are each 
read out by a single PMT. The counters are held in protective aluminum containers with 
a steel fastener, on one corner, which mates to the magnetic shield of the PMT. A few of 
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the counters have special space constraints and are made with wavelength shifting fiber so 
as to allow more flexibility in the orientation of the PMT and magnetic shield. The same 
kind of LED pulser calibration system used in the central muon scintillators is used for the 
Pixel counters. 

The counters are assembled into octants of ~ 100 counters each. The octants provide 
mechanical support for the counters and their infrastructure. These are mounted directly 
onto the A- and B-layer sides of the forward toroid magnets and onto the inside face of the 
C-layer support frame. Two octants are shown in Fig. 18.101 

Forward muon triggers are formed from coincidences in the three layers of scintillation 
counters consistent with a muon of a given momentum. The MDT's provide a pt measure- 
ment at trigger level, especially important in the fiducial region where the fiber tracker has 
reduced coverage. 

8.6.5 Shielding 

The purpose of shielding is to shield the muon detectors from backgrounds generated at 
high | r/\ from the interaction of the beam jets with forward elements of the detector and 
accelerator hardware, such as the beampipe and low-beta quadrupole magnet. The shielding 
is built in several large, movable sections extending from the endcap calorimeters, through 
the forward toroid magnets, to the Tevatron tunnels. The shields themselves totally contain 
the accelerator elements within the collision hall, including the low-beta quadrupole magnet, 
inside a case of 20 inches of iron, six inches of polyethylene, and two inches of lead. The 
shielding is shown in Figs. 18.51 18.61 

8.7 Trigger 

One of the defining features of hadron collider physics is the necessity to select the small 
fraction of all bunch crossings containing interesting collisions. A fast selection process, 



8.7. Trigger 



88 



called the "trigger", sorts events into categories of various levels of interest. An event that 
passes the trigger is written to magnetic tape for later analysis. An event that fails the 
trigger is lost forever. An event can fail to pass the trigger because the collision was a well 
understood process, because the event was mistaken for a well understood process (trig- 
ger inefficiency), or because the trigger or data acquisition were busy processing previous 
collisions (dead time). 

There are three trigger levels. After Level- 1 trigger, the event rate is reduced from 7 
MHz to 10 KHz. Level-2 trigger reduces the event rate 10 times more. We can further 
reduce the event rate 20 times by the Level-3 software trigger. After Level-3, 50 events are 
written to tape per second. 

8.7.1 Trigger Level-1 

The first element of the trigger is the formation of trigger primitives. These are collections 
of a few bits of data that represent the status of various detector elements. For example, one 
set of trigger primitives are the calorimeter cells with energy above a particular threshold. 
Another would be which muon chambers have detected a particle, and whether this particle 
is consistent with a high pt muon, a low pt muon, or neither. Another would be the number 
and position of tracks found in the outer tracker (silicon is not in the trigger at this stage). 

A limited amount of processing is then applied to these trigger primitives. Typically, 
this is performed in Field Programmable Gate Arrays (FPGA's) with inputs from the 
front-end detector electronics. In D0, for instance, Level-1 muon primitives are formed 
from combinations of in-time hits in scintillation counters, coincidences of hit cells in the 
drift chambers, and tracks formed from hit patterns in the axial scintillating fibers of the 
central tracker. 

The primitives are sent to the global Level-1 trigger. Combinations of trigger primitives 
are compared against a runtime programmable list containing the definitions of triggers to 
be used in the run. The D0 Level-1 trigger system can support up to 128 different unique 
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triggers. If the trigger primitives satisfy at least one of the triggers, the event is passed to 
the next trigger level. A block diagram of D0's trigger is shown in Fig. 18.111 



Decision times: -4.2 \ls —100 \ls 

-Detector— LI Trigger— L2 Trigger 
2.3 (7.S) MHz S kHz 



LUM 



FPS 



CFT/ 
CPS 



SMT 



LI 

LUM 



CAL 




L1CAL 


* 



L2Cal — i 



LI 
FPS 



— 



LI 
CTP 



L1CTT 



L2PS 



L2CFT 



L2STT - 



Muon 



LI 

Muon 



L2 
Muon 



—SO ins 

J 1 — L3 Trigger — J 

1 kHz SO Hz 



FPD 




L1FPD 


^ 



L2FW:Combined 
objects (e, jx, j) 



Global 
L2 




L3/DAQ 




Tape 



L1FW: towers, tracks, correlations 



Figure 8.11: Design of the D0 trigger system. 

It takes time to form the primitives and to make a Level-1 trigger decision. Meanwhile 
collisions continue to occur. Since the trigger needs to be ready for the next crossing while 
it is processing an event, the data from the detectors is placed in a pipeline (a micropro- 
cessor's "assembly line" for executing program instructions; a pipelined function unit in a 
processor separates the execution of an instruction into multiple stages) which increments 
every crossing. 

The D0 Level-1 trigger has a deadtimeless output rate of 10 kHz. If any single one of 
the 128 combinations is positive, and the DAQ system is ready for acquisition, then the 
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Level- 1 framework issues an accept, and the event data is digitized and moved into a scries 
of 16 event buffers to await a Level-2 trigger decision. 

The central Level-1 trigger logic is performed locally in the detector octants. A muon 
low pt trigger is defined using only centroids found in the A-layer, while a high px trigger is 
defined by using correlations between centroids found in the A-layer and B- or C-layer. Four 
thresholds (2, 4, 7 and 11 GeV/c) are defined using the CFT information. The information 
for each octant in each region is combined in the muon trigger manager, which produces 
global muon trigger information. The muon trigger manager makes a trigger decision based 
on the pt threshold (2, 4, 7 and 11 GeV/c), pseudorapidity region (\rj\ < 1.0, < 1.5 and 
\rj\ < 2.0), quality (Loose, Medium and Tight) and multiplicity information. This trigger 
decision is sent to the Level-1 Trigger Framework where it is included in the global physics 
trigger decision. 

8.7.2 Trigger Level-2 

At Level-2, the Level-1 decision is confirmed, or not, using the additional time to provide 
more and better information. This is done in two stages. The first stage, called the prepro- 
cessor stage, uses FPGA's to refine the trigger primitives found by Level-1 and to prepare 
the data for the second stage. The second stage, called the global processor stage, combines 
and correlates information from the trigger primitives. 

The preprocessor FPGA's have about 50 /is to perform refinements to the trigger prim- 
itives. For example, the first stage photon primitive, formed from a single tower at Level-1, 
is required to have a shape that looks like an isolated electromagnetic shower at Level-2. 
The additional time is used to compare the energy of the tower of interest with nearby 
towers. Additional related information can also be used, such as energy in the shower max- 
imum. Because of the increased accuracy of the measurement, thresholds can be set tighter 
in Level-2, producing a great deal of the rejection in the global processor stage. 

The global processor stage for D0 is a VME-based CPU card using a Digital Alpha 
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processor. Using a general purpose CPU at this level of the trigger provides a great deal 
of flexibility - essentially all of the Level- 1 primitives are available at Level-2, where they 
can be combined in ways not possible at Level- 1. For instance, Level- 1 can produce photon 
triggers and it can produce jet triggers. However, Level-2 provides the ability to correlate 
the two in ways not possible by the dedicated hardware of Level-1. For example, requiring 
that the jet be opposite the photon in (p and having similar Et- This can be used to select 
events of a given topology with particularly interesting kinematics, for example. Perhaps 
more importantly, though, is that this design provides added flexibility for implementing 
new triggers during the Run II even if they have not been anticipated before the run begins. 
This might be in response to better than promised luminosity from the accelerator, or it 
might be in response to early physics results. 

Level-2 is where the information from the silicon vertex trigger (SVT) becomes avail- 
able. Layers of silicon are read out into this trigger, which uses an associative memory and 
roads provided from the outer tracker to identify silicon tracks. This provides improved 
momentum resolution, but more importantly, it also provides impact parameter resolution 
for each track sufficient to identify particles from displaced vertices. Since the vast majority 
of two track triggers are not from heavy flavor decay, the SVT provides three orders of mag- 
nitude rejection. Silicon information doesn't just improve the impact parameter resolution, 
it also improves the momentum resolution, because the two tend to be highly correlated. 
The maximum deadtimeless output rate of Level-2 is about 1000 Hz at D0. 

By shifting a 3-tube wide window over all the cells in an octant, and looking for wire 
triplets with a matching scintillator hit, the muon preprocessor first finds track stubs sepa- 
rately in the A-layer and the BC-layer. The track stubs found are reported to an ALPHA 
preprocessor board that matches track stubs in the A layer with that in the BC-layer, and 
creates Level-2 objects from matched or unmatched stubs. These Level-2 objects hold the 
(j), i] and pt of the muon, and are reported to the Level-2 global processor. Upon a Level-2 
Accept, the Level-2 objects are sent to Level-3 for more precise muon track reconstruction. 
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8.7.3 Trigger Level-3 

The third level trigger is often described as an event filter. It is a software-based system 
characterized by parallel data paths that transfer data from the detector front-end crates 
to a farm of processors. It reduces the input rate of 1 kHz to an output rate of 50 Hz. 

D0's farm of 500 parallel commodity CPU's builds the event into the offline format, 
runs a modified (for faster execution time) version of the offline reconstruction on the event, 
and makes a decision to accept or reject the event. If the event is accepted, it is already 
in or close to the proper format for offline analysis. Additionally, this trigger level can 
be used to characterize the event and decide whether an event should receive priority in 
reconstruction. Reconstructing a small fraction (say 1%) of events in an "express stream" 
can be used to provide rapid feedback on the detector's performance and . The overall 
output rate of Level-3 is about 50 Hz, with some variation depending on luminosity and 
dataset selection requirements. 

Running what is effectively the offline reconstruction online also provides an excellent 
monitor of the health of the experiment. The full offline reconstruction lags a day or two 
behind the data taking in order to use the final calibration constants, but the Level-3 
reconstruction lags only a fraction of a second. Serious problems can therefore be detected 
before a large amount of data is collected. As an additional benefit, because Level-3 looks 
at the output of Level-2 and offline looks at the output of Level-3, monitoring at Level-3 
examines many times more events than offline monitoring, also improving the probability 
to spot trouble sooner. 

Using commodity processors has a number of advantages. First, the nature of an event 
filter naturally lends itself to parallel processing: each CPU processes an entire event, with 
a supervisor process assigning incoming events to CPUs that finish their events and become 
ready for new ones. Second, these computers are inexpensive, and getting more so with 
time, and finally, the system is highly expandable: additional CPU's can be added at a 
later date, until the bandwidth into Level-3 becomes the limiting factor. 
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The muon Level-3 trigger utilize some aspects of the offline muon reconstruction. Level- 
3 muons have more complete information on the vertex and inner tracking components 
that will improve momentum resolution, and has the ability to require that multiple muons 
come from the same vertex. Similar fits are done in Level-1 as in the final offline recon- 
struction. Requirements on matching the muon track to the inner tracking can reduce 
remnant combinatorics plus punchthroughs. Level-3 also uses the calorimeter information 
to reduce combinatorics, and separate muons into isolated and non-isolated. Level-3 im- 
proves on Level-2's ability to separate muon sources into prompt, slow, or out-of-time by 
fitting the available scintillator hits along a track to the particle's velocity. Level-3 can 
remove a cosmic ray muons both by their being out-of-time and by looking for evidence of a 
penetrating track on the opposite side of the detector. Level-3 can also clean up single muon 
events that Level-1 and Level-2 identified as dimuons, such as those which pass through the 
FAMUS-WAMUS overlap region. 

8.7.4 B Triggers with D0 

S-hadron observability depends strongly on the detector capabilities to trigger on soft 
lepton(s) present in semi-leptonic channel or on J/ip 7 s produced in B decays. Hadronic B 
triggers are not considered in the following. 

The Level-1 muon hardware trigger is based on the combination of low px track candi- 
dates measured in the CFT, spatially matched with hits in the scintillator planes and/or 
drift chambers. Single muon events with p T > 4 GeV/c run prescaled currently, and di- 
muons with p T > 2 GeV/c are expected to run unprescaled. 

The electron trigger is aimed at soft electron pair detection. Level-1 candidates are 
selected separately in EM calorimeter trigger towers (Arj = A(p = 0.2) with a transverse 
energy deposit Et > 2.0 GeV, and in the tracking system with a low px track coincident 
with pre-shower cluster. Electron candidates of both systems are then required to match 
within a quadrant in ip and to have opposite signs. These electron triggers are not used in 
the following analysis. 
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Level-2 triggers include setting of the invariant mass window and angular cuts in di- 
lepton channels to select J ftp decays and improve background rejection. 

8.8 Data Acquisition 

The data flow in trigger Level-3/data acquisition is illustrated in Fig. 18.121 Events are 
reconstructed on the FNAL processor farm system, with that portion dedicated to D0 
capable of matching the 50 Hz data acquisition rate. Following reconstruction, data are 
stored on a tightly coupled disk and robotic tape system, and made available for analysis 
on a centralized analysis processor. 
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Figure 8.12: Simplified data flow in Level-3/DAQ for D0 Run II. 



Chapter 9 



Establishing the A£ Signal 

9.1 Data and Monte Carlo Event Samples 

The data used corresponds to an integrated luminosity of approximately 225 pb" 1 col- 
lected between April 2002 and January 2004. This analysis is described in Ref. [66] and 
repeated in Ref. [67]. 

9.1.1 Dimuon data sample 

The data used in this analysis was part of a skim that was called the "D" skim, which 
included either D mesons for semileptonic studies or J/^'s for exclusive channels. The 
dimuon requirements on this skim are: 

• two muons with opposite charge; 

• pt > 1 GeV for each muon; 

• number of CFT hits > 1 for each particle; 

• at least 1 muon with nseg = 3; nseg is a muon quality cut that describes the hits in 
the muon detector layers. 

• for nseg > 0, Pt{h) > 1-5 GeV; for nseg = 0, pt > 1 GeV; 

• 2.5 GeV < MassOiV") < 3.6 GeV; 
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• for a muon with nseg = 0: 

- pt(/x + /x") > 4.0 GeV; 

- p(ji) < 7 GeV; 

- pr of second muon > 2.5 GeV; and 

- x 2 °f global muon track fit < 25 for both muons. 

9.1.2 Primary vertex reconstruction 

To determine the primary vertex, the AATrack package [59] is used: 

AATrack vO-10-06-01 

beam spot file version 2.09 

that follows the procedure as outlined in Ref. [60] . A primary vertex is determined for each 
event by minimizing a x 2 function that depends on all the tracks in the event and a term 
that represents the beam spot constraint. The beam spot is the run-by-run average beam 
position, where a run typically lasts several hours. The beam spot is stable during the 
periods of time when the proton and the antiproton beams are kept colliding continously, 
and can be used as a constraint for the primary vertex fit. The initial primary vertex 
candidate and its x 2 are obtained using all tracks. Next, each track used in the x 2 calculation 
is removed temporarily and the x 2 is calculated again; if the x 2 decreases by 9 or more, this 
track is discarded from the PV fit. This procedure is repeated until no more tracks can be 
discarded. The PV resolution achieved by using this procedure is 25 fim. 

9.1.3 Secondary vertex reconstruction 

To reconstruct secondary vertices (i.e., J ftp, A and A° vertices), the BANA package [61] is 
used. 
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9.1.4 Ag Candidate Selection 

Further cuts are used to select J /if) mesons: 

• Require that a primary vertex reconstructed. 

• For dimuons forming a vertex, require: 

— distance from primary vertex to J /if) vertex < 10 cm. 

— Limited number of track measurements downstream (< 1) of the vertex and 
missed hits upstream (< 5) of the vertex. 

• 2.8 GeV < Mass(/i+/x~) < 3.35 GeV; and 

• number of SMT hits for each muon >= 1. 

Figure 19.11 and Fig. 19.21 show the fj,~) invariant mass for all inclusive events and for 
events containing a A candidate, respectively. 

To select A — ► pir~ candidates, we require: 

• two oppositely charged tracks forming a vertex; 

• the higher momentum track is assumed to be the proton; 

• limited number (< 1) of track measurements downstream of the vertex; 

• number of CFT hits > 1 for each track; 

• 1.105 GeV < Mass(pvr) < 1.125 GeV; 

• p T (A) > 0.4 GeV; and 

• A candidates falling within a Kg mass window of 0.465 GeV < Mass(7r + 7r~) < 
0.52 GeV are removed. 
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Figure l9~3l shows the invariant mass of (pir) for all events containing a J ftp candidate. 

Once good J/ip and A candidates have been found in the same event, the two muons 
from J ftp are combined with the reconstructed neutral A track to form a Ajj vertex. Further 
requirements are made on the vertex: 

• vertex x 2 (^b) < 10; 

• p{A° b ) > 5 GeV; 

• xy distance from J/ip vertex to A vertex > 0.3 cm; 

• assuming that the higher momentum track is the proton, p{p) > 0.8 GeV; and 

• collinearity of A® > 0.99, where collinearity is defined as the cosine of the angle between 
the momentum of A^ in the xy plane and the direction from primary to secondary in 
the xy plane. 

Figure l9~H shows the invariant mass of (/U + /i~A) subject to the above cuts. The signal is 
modeled with a Gaussian function, and the background by a constant plus an exponential 
function. The exponential function is taken for the fact that there is a large shoulder on 
the left of the signal due to partially reconstruced B mesons. The number of signal events 
extracted from the fit is 48±11 with a mean of 5630±10 MeV and a width of a = 39±8 MeV. 

To check the validity of the signal for the presence of long lived signal for lifetime 
measurement purposes, the decay length cut and decay length significance cut were varied. 
Figure l9~51 shows the decay length and decay length significance in data. 

A transverse signed decay length is defined using the vector between the A® vertex and 
the primary vertex, L xy = x Ab - x PV , and forming: L xy = (L xy ■ p T (A&))/|pr(A 6 )|. 

Figures 19.61 19.71 and 19.81 show the invariant mass of the /i + /x~A system subject to the 
above cuts and adding L xy > 0.01 cm, L xy > 0.02 cm, and L xy > 0.03 cm requirements, 
respectively. The number of signal events extracted from the fit in each case is 33 ± 8, 
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33 ± 7, and 31 ± 6, respectively, with the signal-to-noise increase in each case as expected 
for a signal with lifetime. 

Figures 19. 91 and f9 . 101 show the invariant mass of the /x + /i~A system subject to the above 
cuts and adding L xy / a{L xy ) > 2 and > 4, respectively. The number of signal events 
extracted from each fit is 37 ± 8 and 31 ± 6, respectively, with increased signal-to-noise, 
again indicating lifetime information in the signal. 

Figures [9. 1 11 19. 121 and 19 .131 show the invariant mass of the /i + /i~A system subject to the 
above cuts and adding proper lifetime cuts cr > 370 fim, ct > 740 /iin, and cr > 1000 /im 
requirements, respectively. The number of signal events extracted from the fit in each case 
is 18 ± 5, 7 ± 3, and too small to allow a fit, respectively. This fractional event loss is 
consistent with the signal having a lifetime of the order of a B hadron. 
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Figure 9.1: Invariant mass distribution of the system for all the events. The signal is 

described by a double Gaussian function and the background by a second order polynomial. 
The signal mean is 3070.0 ± 0.1 MeV and the fitted widths are a\ = 55 ± 1 MeV and 
02 = 100 ± 1 MeV. 
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Figure 9.2: Invariant mass distribution of the system for all the events containing 

a A candidate. The signal is described by a double Gaussian function and the background 
by a second order polynomial. The signal mean is 3069 ± 1 MeV and the fitted widths from 
double Gaussian function are a\ = 52 ± 1 MeV and 02 = 91 ± 5 MeV. 
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Figure 9.3: Invariant mass distribution of the (p, tt) system for all the events containing a 
J/ifj candidate. The signal is described by a double Gaussian function and the background 
by a second order polynomial. The signal mean is 1115.3 ± 0.1 MeV and the fitted widths 
are a\ = 4.6 ±0.1 MeV and 02 = 2.0 ±0.1 MeV for a double Gaussian function. 
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Figure 9.4: Invariant mass distribution of the (J/ifi, A) system for all A° candidates. The 
signal is described by a Gaussian function and the background by a constant plus an expo- 
nential function. 
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Figure 9.5: Top: decay length distribution; and bottom: decay length significance distribu- 
tion for data. 
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Figure 9.6: Invariant mass distribution of the (J ftp, A) system for all A° candidates plus 
the requirement L xy > 0.01 cm. The signal is described by a Gaussian function and the 
background by a constant plus an exponential function. 
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Figure 9.7: Invariant mass distribution of the (J ftp, A) system for all candidates plus 
the requirement L xy > 0.02 cm. The signal is described by a Gaussian function and the 
background by a constant plus an exponential function. 
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Figure 9.8: Invariant mass distribution of the (J/ip,A) system for all A^ candidates plus 
the requirement L xy > 0.03 cm. ehe signal is described by a Gaussian function and the 
background by a constant plus an exponential function. 
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Figure 9.9: Invariant mass distribution of the (J/ip,A) system for all candidates plus 
the requirement L xy /a(L xy ) > 2. The signal is described by a Gaussian function and the 
background by a constant plus an exponential function. 
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Figure 9.10: Invariant mass distribution of the («//?/>, A) system for all A® candidates plus 
the requirement L xy /a(L xy ) > 4. The signal is described by a Gaussian function and the 
background by a constant plus an exponential function. 
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Figure 9.11: Invariant mass distribution of the (J/ip,A) system for all A® candidates plus 
the requirement cr > 370 fim. The signal is described by a Gaussian function and the 
background by a constant plus an exponential function. 
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Figure 9.12: Invariant mass distribution of the (J/ip,A) system for all candidates plus 
the requirement cr > 740 fim. The signal is described by a Gaussian function and the 
background by a constant plus an exponential function. 
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Figure 9.13: Invariant mass distribution of the (J/ip,A) system for all candidates plus 
the requirement cr > 1000 /an. The signal is described by a Gaussian function and the 
background by a constant plus an exponential function. 
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9.1.5 Monte Carlo event samples 

9.1.6 MC Signal A° -> J/i/jA 

To simulate the decay chain A® — > J/tpA, J/ip A — > £>7r~, we use the evtgen decay 

model and the pythia generation program [63]. Below we show the decay file used. 

j 

DECAY LAMB 
CHANNEL 1.000 PSI 
ENDDECAY 

J 

DECAY PSI 

ANGULAR_HELICITY -1 
ANGUL AR_HEL I C I T Y 
ANGUL AR_HEL I C I T Y 1 
CHANNEL 40 1.000 MU+ 
ENDDECAY 

I 

DECAY LAM 
CHANNEL 1.000 P+ 
ENDDECAY 
> 

Before passing the generated events through the suite of programs for the detector 
simulation, hit simulation, trigger simulation, and track and particle reconstruction, we 
apply the following "pre-GEANT" selection cuts using the dOmess package [64] : 

• presence of the decay chain A° — > J/tpA. 

• pr(n) > 1.5 GeV and 0.8 < |r/| < 2.0 or p T {p) > 3 GeV and 0.8 < \rj\. 



LAM 



1. 0. 
1. 0. 
1. 0. 
MU- 



PI- 
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"GEANT" is a detector simulation tool, and is described in detail in Ref. [65]. 

We apply the kinematic and quality cuts as described in the previous section. After 
kinematic cuts, we end up with 50000 events. Since the A is a long-lived particle, it decays 
at the GEANT processing stage, so out of 50000 events, 31500 events contain the decay 
A — ► pir~ . The same analysis cuts as data are applied. The number of events passing 
all event selection criteria is 315. The reconstruction efficiency after the kinematic cuts 
is e = (1.0 ± 0.06)%. Figure 15.141 shows the invariant mass of the combination (A, J/VO 
for the Monte Carlo sample. The mass and the width observed in MC are M(A, J/VO 
=5642 ± 2 MeV and a = 35 ± 2 MeV respectively. The mass is slightly lower in data 
compared to MC due to uncertainty in D0 momentum scale, and the fitted width observed 
in data is consistent with that obtained in the MC. 

Figure 15.151 shows the decay length and decay length significance distribution from the 
MC within the mass window of A^. These distributions clearly indicate that the signal 
candidates have a long lifetime component, as expected. 

9.2 B Q d Event Selections 

The decay — > J/tpK® has similar topology to A® — > J/ipA decay and has higher 
statistics, making it a good test sample. Both decays are reconstructed the same way, the 
only difference is Kg is that a reconstructed K s — > 7t + tt~ is combined with the J/tp instead 
of a A. To select K s — > ir + ir~ candidates, we require: 

• two oppositely charged tracks forming a vertex; 

• limited number of track measurements downstream (< 1) of the vertex; 

• number of CFT hits > 1 for each track; 

• 0.465 GeV < Mass(7r+vr-) < 0.52 GeV; 
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• prills) > 0.4 GeV. 

Figure ID. 161 shows the invariant mass of (tttt) for all events containing a J/tp candidate. 

Once good J ftp and K% candidates have been found in the same event, the two muons 
from J/ip are combined with the reconstructed neutral Kg track to form a B® vertex. 
Further requirements are made on the vertex: 

• vertex x 2 (Bd) < 25; 

• p{B° d ) > 5 GeV; 

• collinearity of B^ > 0.99, where collinearity is defined as the cosine of the angle 
between the momentum of B® in the xy plane and the direction from primary to 
secondary vertex in the xy plane. 

Figure ID. 171 shows the invariant mass of n + fj,~K s subject to the above cuts. The signal 
is modeled with a Gaussian function, and the background by a second-order polynomial. 
The number of signal events extracted from the fit is 300 ± 39 with a mean of 5270 ± 5 MeV 
and a width of a = 38 ± 5 MeV. 

Figures 19. 18(19.191 and 19.201 show the invariant mass of the fi + , K s system subject to 
the above cuts and adding L xy > 0.01 cm, L xy > 0.02 cm, and L xy > 0.03 cm requirements, 
respectively. The number of signal events extracted from the fit in each case is 299 ± 30, 
298 ± 25, and 278 ± 22, respectively. 

Figures 19.211 and 19.221 show the invariant mass of the fj, + ,fj,~,Kg system subject to the 
above cuts and adding L xy /a(L xy ) > 2 and > 4, respectively. The number of signal events 
extracted from each fit is 300± 27 and 259 ± 21, respectively. 

Figures 19.231 19.24*1 and 19.251 show the invariant mass of the n + ,fi~,Kg system subject 
to the above cuts and adding cr > 460 fim, ct > 920 /mi, and cr > 1320 fim requirements, 
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respectively. The number of signal events extracted from the fit in each case is 153 ± 16, 
60 ± 10, and 22 ± 7, respectively. These numbers are consistent with an exponentially falling 
B lifetime. 
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Figure 9.14: Invariant mass distribution of candidates in MC simulated events. A double 
Gaussian function was used to model the signal, and a constant plus an exponential function 
was used for the background. 
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Figure 9.15: Top: decay length distribution; bottom: decay length significance distribution 
for MC simulated events. 
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Figure 9.16: Invariant mass distribution of the (7r + 7r _ ) system for all the events with a J/tj) 
candidate. The signal is described by a double Gaussian function and the background by 
a second order polynomial. The signal mean is 493.5 db 0.1 MeV and the fitted widths are 
ffi = 5.9 ± 0.1 MeV and a 2 = 11.7 ± 0.1 MeV. 
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Figure 9.17: Invariant mass distribution of the (J/ip, Kg) system for all B^ candidates. The 
signal is described by a Gaussian function and the background by a second order polynomial. 
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Figure 9.18: Invariant mass distribution of the (J/ip^ Ks) system for all candidates plus 
the requirement L xy > 0.01 cm. The signal is described by a Gaussian function and the 
background by a second order polynomial. 
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Figure 9.19: Invariant mass distribution of the (J/tp, Kg) system for all candidates plus 
the requirement L xy > 0.02 cm. The signal is described by a Gaussian function and the 
background by a second order polynomial. 
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Figure 9.20: Invariant mass distribution of the (J/ifi, Ks) system for all candidates plus 
the requirement L xy > 0.03 cm. The signal is described by a Gaussian function and the 
background by a second order polynomial. 
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Figure 9.21: Invariant mass distribution of the (J/tp, Kg) system for all candidates plus 
the requirement L xy /a(L xy ) > 2. The signal is described by a Gaussian function and the 
background by a second order polynomial. 
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Figure 9.22: Invariant mass distribution of the (J/ifi, Kg) system for all B® candidates plus 
the requirement L xy /a(L xy ) > 4. The signal is described by a Gaussian function and the 
background by a second order polynomial. 
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Figure 9.23: Invariant mass distribution of the (J/tp, K s ) system for all B^ candidates plus 
the requirement cr > 460 ^m. The signal is described by a Gaussian function and the 
background by a second order polynomial. 
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Figure 9.24: Invariant mass distribution of the (J/ifi, Ks) system for all candidates plus 
the requirement cr > 920 /im. The signal is described by a Gaussian function and the 
background by a second order polynomial. 
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Figure 9.25: Invariant mass distribution of the (J/tp, Kg) system for all candidates plus 
the requirement cr > 1320 ^m. The signal is described by a Gaussian function and the 
background by a second order polynomial. 
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Aj) and Bd Lifetime Measurement 

10.1 Proper Decay Length 

We define the signed decay length of a baryon as the vector pointing from the primary 
vertex to the decay vertex projected on the A° momentum in the transverse plane: 

A 

Lxy = {x A o - x prim ) -pr/PT, (10.1) 

where pr is the measured transverse momentum vector and pt is its magnitude. The 
primary vertex is reconstructed individually for each event as described earlier. 

The proper lifetime, r, and the proper decay length, cr, are then defined by the relation: 

A 

ct = L x *-M a o/pt, (10.2) 

b 

where M^p = 5624 ± 9 MeV is the world average mass of the hP h baryon [29]. 



10.1.1 Proper Decay Time Distributions 

We divide the A° mass range into three bands: "low side", 5.37 < M(A°) < 5.5 GeV; 
"middle", 5.5 < M(Ag) < 5.73 GeV; and "high side", 5.73 < M(A°) < 5.87 GeV. The 
middle mass band includes the signal; all three contain comparable numbers of background 
events. The proper lifetime distributions in the three mass bands are compared in Fig. 111). fl 
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The distributions are normalized to the number of events in the signal region. The low and 
high sidebands are dominated by physics background, such as B meson background. We 
assume the same parametrization of the background shape in the entire mass region. 
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Figure 10.1: The proper lifetime distribution in three mass bands (see text). 



10.2 Fitting Procedure for A° b 



We use the MINUIT minimization program in root (TMinuit class). The source code 
and user tips are available at [69] . 

The proper decay length and invariant mass distribution of the A° candidates are fit 
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simultaneously using an unbinned maximum log-likelihood method. The likelihood function 
C is given by: 



where iV=353 is the total number of events, J- l sig is the product of the signal mass and proper 
decay- length probability density functions, Fl ck is the product of the background mass and 
proper decay-length probability density functions, and f S i g is the fraction of signal. 

We use a range of 5.0 GeV < M(A^) < 6.3 GeV for the mass window. Previously we 
implemented lifetime related cuts (collinearity and decay length related cuts), which were 
used to solidly establish the signal. However, for lifetime purposes we have to remove these 
cuts. We will use the signal with no lifetime cuts as shown in Fig 111). 2l The cuts used for 
selecting these events are: 

• vertex x 2 (M) < 12; 



• xy distance from J ftp vertex to A vertex > 0.9 cm; 

• p T (A) > 2.4 GeV; and 

• cr(cr) < 100 fs. 

These are the events used to extract the lifetime. Figure I1U.2I shows the invariant mass 
of fi + fi~A subject to the above cuts. The signal is modeled with a Gaussian function, 
and the background by a constant plus a second order polynomial. The number of signal 
events extracted from the fit is 32 ± 10 with a mean of 5622 ± 12 MeV and a width of 
a = 36 ± 10 MeV. 

The lifetime distribution of the signal is parameterized by an exponential convoluted 
with a Gaussian function. The lifetime resolution of background is approximated by a 




(10.3) 



• p(Ag) > 5 GeV; 
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("^pntral value -\- error 


Proper decay length, cr 


4Uo ± 59/jm 


Signal fraction, f s %g 


u.uyi it U.Uzo 


Mass of A 6 , M(A b ) 


5622±12 MeV 


Width of Aft, o-(Aft) 


36± 10 MeV 


Mass Slope(O), aO 


-0.31 ±0.031 


Mass Slope(l), al 


0.026 ± 0.028 


Error scale factor, s 


1.7 ± 0.1 


slope — bkg — neg 


88 ± 52 /im 


slope — bkg — pos 


423 ± 114 fim 


norm — bkg — neg 


0.02 ± 0.02 


norm — bkg — pos 


0.1 ± 0.02 



Table 10.1: Unbinned maximum likelihood fitting results for the A^ signal. 



superposition of a Gaussian function centered at zero, one exponential for the negative cr 
region and one exponential for the positive cr region, with free slopes. The width of the 
Gaussian functions is taken from the event-by-event measurement of decay length error, o~l. 
To allow for the possibility of the lifetime uncertainty to be systematically mis-estimated, 
we introduce a free scale factor s that multiplies the decay length error, i.e., s ■ o~l. The 
fraction of the signal is set to the central value obtained in the separate mass fit and allowed 
to vary within its statistical errors. This fraction is f s i g = N s i gna i/N total in the defined mass 
window. 



Fit to the data 



The fit results are summarized in Table fTTTTTl 



The lifetime distribution with the fit results overlaid is shown in Fig. 110.31 In the figure, 
for illustration purposes, we assume a constant decay length resolution of 50 /im, which is 
the average value observed in the data. 

In Fig. 110.41 we fit only the sidebands with the background functional form showing a 
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good fit as expected. 



For the mean proper decay length we obtain: 



ct(A") = 408 ± 89 (stat) /im. 



(10.4) 



The central value of mean proper decay length is larger than what has been reported 



previously [29] (367 ± 24 /im); however, it is consistent within statistical errors. 

The fraction of signal in the sample is f s i g = 0.091 ± 0.028. The fitted value of the 
lifetime resolution scale factor, s = 1.7 ± 0.1, is inconsistent with unity which means that 
our errors in the data have been underestimated. This has been shown also in other lifetime 
analyses [36]. 

10.3 Fitting Procedure for Bd 

We used the same fitting procedure for Bd as for A°. The proper decay length and 
invariant mass distribution of the candidates are fit simultaneously using an unbinned 
maximum log-likelihood method. The likelihood function C is given by: 



where A r =2409 is the total number of events, T\ ig is the product of the signal mass and 
proper decay-length probability density functions, ^ cfc is the product of the background 
mass and proper decay-length probability density functions, and f s i g is the fraction of signal. 

We use a range of 4.9 GeV < M(-B<f) < 5.6 GeV for the B^ mass window. We will use 
the signal with no lifetime cuts as shown in Fig 110.51 The cuts used for selecting these 
events are: 




(10.5) 



• vertex x 2 (Bd) < 25; 



10.3. Fitting Procedure for B d 



128 





(Central value -\- error 


Proper decay length, ct 


4J8 ±oa (im 


Signal fraction, f s ig 


U.loO It U.U10 


Mass of B d , M(B d ) 


5273 ± 6 MeV 


Width of B d , o{B d ) 


50 ± 6 MeV 


Mass Slope (0), aO 


-0.37 ±0.003 


Mass Slope (1), al 


0.034 ± 0.003 


Error scale factor, s 


1.56 ± 0.05 


slope — bkg — neg 


131 ± 19 fim 


slope — bkg — pos 


466 ± 35 fim 


norm — bkg — neg 


0.040 ± 0.008 


norm — bkg — pos 


0.18 ± 0.01 



Table 10.2: Unbinned maximum likelihood fitting results for the B^ signal. 

• p(B d ) > 5 GeV; 

• xy distance from J/i/j vertex to K s vertex > 0.3 cm; and 

• p T {K%) > 1.8 GeV. 

These are the events used to extract the lifetime. Figure 110.51 shows the invariant mass 
of (// + /U _ A°) subject to the above cuts. The signal is modeled with a Gaussian function, 
and the background by a constant plus a second order polynomial. The number of signal 
events extracted from the fit is 330 ± 39 with a mean of 5622 ± 12 MeV and a width of 
a = 50 ± 6 MeV. 

Signal and background are parametrized identically to A^. 

Fit to the data for B d 

The fit results are summarized in Table HTPl 



The lifetime distribution with the fit results overlaid is shown in Fig. 110.61 In the figure 
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for illustration purposes, we again assume a constant decay length resolution of 50 /im, 
which is the average value observed in the data. 

For the B d mean proper decay length we obtain: 

cr{B d ) = 428 ± 35 (stat) (jm. (10.6) 

The central value of B d mean proper decay length is in a good agreement with previous 
measurements of this value [29] and the world average of cr(B d ) = 460 ± 4 (stat) //m. 

The fraction of signal in the sample is f s i g = 0.135 ± 0.016. The fitted value of the 
lifetime resolution scale factor, s = 1.56 ± 0.05, is consistent with that observed for the A? 
channel. 

10.4 Systematic Errors for A5 and Bd Lifetime 

We summarize the systematic errors for our measurements for A|j and B d lifetime in 
Table ITTH1 

The contribution from the uncertainty in detector alignment is estimated by reconstruct- 
ing the B° sample with the positions of the SMT sensors shifted outwards radially by the 
alignment error of 30 /im in the radial position of the sensors and then re-fitting for the 
lifetime. 

To check the stability of the fits, we varied the fitted parameters other than cr within 
±<t of their nominal value, one at a time. This systematic turned out to be the dominant 
error for A° but, comparable with the other ones for B d . 

To estimate the uncertainity in the modeling, we used different background shapes to de- 
scribe the background. We used one and two exponential PDF for the lifetime background. 
We also used an exponential and second order polynomial PDF for the mass. 

We also fitted the background only separately and the fitted parameters from this fit are 
used as fixed parameters when we fit background plus signal. The difference in the results 
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Source 


A 6 (fi m) 


B d (/im) 


Ratio 


Alignment 


5 


5 


0.016 


Variation of Fitting Parameters 


16 


5 


0.039 


Modeling 


11 


3 


0.026 


Fixing the parameters from sideband fit 


1 


9 


0.020 


Total 


20 


12 


0.050 



Table 10.3: Summary of systematic uncertainties in the measurement of cr for A® and B d 
and their ratio. The total uncertainties are also given combining the individual uncertainties 
in quadrature. 



between the lifetime we get when using fixed background and the unfixed background gives 
rise to only a small uncertainity. 
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Figure 10.2: Invariant mass distribution of the (J/ifi, A) system for all A® candidates without 
any lifetime cuts. The signal is described by a Gaussian function and the background by a 
second order polynomial. 
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Figure 10.3: The proper decay length, cr, of the A° candidates. The curves show: the signal 
contribution, dotted (green); the background, dashed (red); and total, solid (black). 
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Figure 10.4: The proper decay length, cr, of the sidebands. 
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Figure 10.5: Invariant mass distribution of the (J/ip,Ks) system for all candidates 
without lifetime cuts. The signal is described by a Gaussian function and the background 
by a second order polynomial. 
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Figure 10.6: The proper decay length, cr, of the candidates. The curves show: the 
signal contribution, dotted (green); the background, dashed (red); and total, solid (black). 



Chapter 11 



Discussion and Conclusion 

11.1 Summary on Establishing the A® Signal 

We observe evidence for the baryon using the exclusive decay A^ — > J/tpA followed 
by J/ip — ► [i + fi~ and A — > jot". The data sample corresponds to an integrated luminosity 
of approximately 225 pb _1 . We find 48±10 A° candidates and determine the mass of A° 
baryon to be 5630 ± 10 (stat) MeV. This result is consistent with previous A° baryon mass 
measurements [29] . The A° yield following lifetime-related cuts also indicates that the signal 
contains a lifetime distribution different from background. The decay B® — > J/ipK^ is used 
as a test sample, so for completenness, we also show the results for B% meson. We find 
300±39 B® candidates and determine the mass of the B® meson to be 5270 ± 5 (stat) MeV. 

11.2 Summary of Preliminary Lifetimes 

We were able to extract the lifetimes of A^ — ► J/t(jA and B® — > J/t/iKg. The lifetime 
of the A° baryon has been measured to be 

cr (A£) = 408 ± 89 (stat) ± 20 (syst) urn. (11.1) 

r(Ag) = 1.36 ± 0.30 (stat) ± 0.07 (syst) ps. (11.2) 

This result is consistent with the world average [29] of r(A°) = 1.229 ± 0.080 ps. 
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The lifetime of the B d baryon has been measured for completeness: 



CT(B d ) = 428 ± 35 (stat) ± 12 (syst) //m. 



(11.3) 



T(B d ) = 1.43 ± 0.12 (stat) ± 0.04 (syst) ps. 



(11.4) 



This result is in agreement with the world average [29] of T{B d ) = 1.536 ± 0.014 ps. 
Finally we report the ratio t(Ai,)/t(B ( i): 

444u P = °- 95 ± °- 22 ( stat ) ± °- 05 ( s y st )' ( n - 5 ) 

which is consistent with previous world measurements and with the theoretical predictions. 
Although the data used in this study are limited we look forward to adding more data and 
should be able to carry out a competitive measurement of the A° lifetime via this exclusive 
decay channel. 
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Further Details on Systematics Error 
Determination 

We will show here a case study for systematics: we assume the mixing parameters x and 
y are 5% and Im £ = 0, and investigate Re £. We previously observed that MC corrections 
were very small in what we call the standard FOCUS absorption model. It includes D° 
and D° cross sections half the pion cross section, Moliere scattering and the absorbed 
particles having no byproducts. While interactions of pions and kaons with matter have 
been measured, no such measurement exist for intereactions of charm particles with matter. 
We assumed that a(D°) = a(D°) = 1/2 x cr(vr) and even though this is a good assumption 
we want to examine the sensitivity of CPT parameters when cross sections are not equal. 
In Fig. lA.lf a) several variations from standard absorption are shown. We observe that we 
are not very sensitive to the difference between the cross sections of D° versus D°. For 
instance, increasing the D° cross section by 20% relative to D°, results in a 4% decrease 
in Re £. In reality these cross sections should be closely equal which results in systematic 
effect of less than 1% associated with typical uncertainties in charm cross sections. 

Fig. lA.lf b) shows the stability of Re £ versus split samples. Some momentum depen- 
dence is observed but the variations are not very significant. 

We vary — and Kaonicity (W n — Wk) requirements. We observe in Fig |A.2f b) the fitted 
Re£ versus 16 combinations of — and (W^ — Wk) for x and y mixing of 5%. We saw that 
the corrections remain small independent of the set of cuts we used. From the way these 
cuts are applied there is a correlation of the data from one set to the other one but the 
underlying systematic effect is small. 

We also investigated bin and sideband variations to check for any systematic effect. 
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Figure A.l: Re £ for different MC absorption models. The order of each point from left 
to right is 1) uncorrected 2) standard FOCUS (baseline): 3) baseline and byproducts 4) 
baseline and charm absorption zero 5) baseline, charm absorption zero and byproducts, 6) 
baseline, a(D°) = 1/2 x a(ir + ) = 11.5mb and a(D°) = l/2_x cr(vr~) = 13m6, 7) baseline 
and a{D°) = 3/2 x a(D°), 8) baseline and a(D°) = 3 x a(D°) (a). Different split samples 
(b). 
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Figure A. 2: Stability of Ret; for (a) different fit variants and (b) cut variants (^,Kaonicity) 
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Fig lA.2r a) shows five different sideband methods and two different bin widths. The five dif- 
ferent sidebands are (±2<t,±3<t,±7ct), (±2ct,±4<t,±8<t), (±2<r,±3cr,±9cr), (±2<7,3<7,7<r) and 
(±2cr,— 3er,— 7a) and two bins are 500 fs and 300 /s. In Fig |A.2b ). the first 5 points corre- 
spond to 500fs bin width and the last five points correspond to 300 fs bin width. 
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